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ABSTRACT 

We construct a simple model of the star-formation- (and resultant supernova-) 
driven mass and energy flows through the inner ~200 pc (in diameter) of the Galaxy. 
Our modelling is constrained, in particular, by the non-thermal radio continuum and 
7-ray signals detected from the region. The modelling points to a current star-formation 
rate of 0.04 — 0.12MQ/year at 2a confidence within the region with best-fit value in 
the range 0.08 — 0.12 M0 /year which - if sustained over 10 Gyr - would fill out the 
~ 10^ Mq stellar population of the nuclear bulge. Mass is being accreted on to the 
Galactic centre (GC) region at a rate AIjn ^ 0.3 Mq /year. The region's star-formation 
activity drives an outflow of plasma, cosmic rays, and entrained, cooler gas. Neither 
the plasma nor the entrained gas reaches the gravitational escape speed, however, 
and all this material fountains back on to the inner Galaxy. The system we model 
can naturally account for the recently-observed ^ 10^ Mq 'halo' of molecular gas 
surrounding the Central Molecular Zone out to 100-200 pc heights. The injection of 
cooler, high-metallicity material into the Galactic halo above the GC may catalyse the 
subsequent cooling and condensation of hot plasma out of this region and explain the 
presence of relatively pristine, nuclear-unprocessed gas in the GC. This process may 
also be an important ingredient in understanding the long-term stability of the GC 
star-formation rate. The plasma outflow from the GC reaches a height of a few kpc 
and is compellingly related to the recently-discovered Fermi Bubbles by a number of 
pieces of evidence. These include that the outflow advects precisely i) the power in 
cosmic rays required to sustain the Bubbles' 7-ray luminosity in saturation; ii) the hot 
gas required to compensate for gas cooling and drop-out from the Bubbles; and iii) 
the magnetic field required to stabilise the walls of these structures. Our modelling 
demonstrates that ~ 10^ Mq of hot gas is processed through the GC over 10 Gyr. 
We speculate that the continual star-formation in the GC over the age of the Milky 
Way has kept the SMBH in a quiescent state thus preventing it from significantly 
heating the coronal gas, allowing for the continual accretion of gas on to the disk and 
the sustenance of star formation on much wider scales in the Galaxy. In general, our 
investigations explicitly reveal the GC's important role in the Milky Way's wider stellar 
ecology. 

Key vifords: cosmic rays - galaxies: star formation - Galaxy: center - ISM: jets and 
outflows - ISM: supernova remnants - radiation mechanisms: non-thermal 



The inner 200 pc (in diameter) of the Milky Way features a 
spectacular confluence of unusual and energetic astrophysical 
phenomena. Within this region of the Galaxy - circumscribed 
by the Inner Lindblad Resonance associated with the non- 
axisymmetric gravitational potential of the Galactic bar 
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1 INTRODUCTION 



- the distribution of stars cusps sharply into the distinct 
population of the so-called nuclear bulge (Serabyn & Morris 
1996). Correspondingly, over the same region the current, 
inferred areal density of star formation, E«, sharply peaks 
to ^ 200 A/o/kpc^/yr. This is approximately three orders of 
magnitude higher than the mean value in the Galactic disk. 
With such a high E*, observations of the nuclei of external, 
star-forming galaxies tell us to expect a star-formation-driven 
outflow; there is much empirical and theoretical evidence 
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that such an outflow exists in the GC as we have explored 
in a number of recent papers (Crocker et al. 2011a; Crocker 
& Aharonian 2011; Crocker et al. 2011b). This paper adds 
significantly to that evidence. 

The very high star formation rate (SFR) density like- 
wise sustains a very high energy density in all phases of the 
GC ISM. Most directly, the optical and UV output of the 
many young, hot stars in the region is reprocessed by thick, 
ambient dust into a dominantly infrared photon background 
of ~100 eV cm""^. Radio continuum and 7-ray observations 
allow one to place a robust lower limit of ~50 pG on the 
typical magnetic field throughout the entire inner ~80Q pc 
(in diameter) of the Galaxy (Crocker et al. 2010a); modelling 
(Crocker et al. 2011a,b) points to a magnetic field in the 
inner 200 pc that is at least 100 fiG. In association with 
and, in fact, as a necessary precondition to, the high SFR, 
observations reveal an enormous agglomeration of hot, dense, 
and highly-turbulent molecular gas of mass 3 X lO'' Mq (Dah- 
men et al. 1998; Molinari et al. 2011), 5-10% of the Milky 
Way's entire H2 allocation. This gas forms an asymmetric 
distribution extended along the plane referred to as the Cen- 
tral Molecular Zone (CMZ: Serabyn & Morris 1996). Recent 
infrared observations of the CMZ by Herschel (Molinari et 
al. 2011) place much of this H2 on a ~ 100 pc-radius ring, 
apparently akin to the star-forming rings observed in the 
nuclei of many face-on galaxies. 

The CMZ H2 is constantly bombarded by an extended, 
hard-spectrum cosmic ray ion population which results in 
a diffuse glow of hard-spectrum, ~TeV 7-rays coextensive 
with the gas (Aharonian et al. 2006). GC X-ray observations 
(Koyama et al. 1989) apparently reveal the existence of a 
very hot (~ 7 keV), extended thermal plasma which would 
have an energy density similar to the lightfield, turbulent 
molecular gas, and magnetic field (Spergel & Blitz 1992). 

It must also be remarked that the GC hosts the Milky 
Way's resident supermassive (~ 4.3 x 10^ Mq: Gillessen et 
al. 2009) black hole (SMBH). Though currently in a state of 
apparently unusual quiescence, this must certainly have been 
much more active at various times in the past (e.g. Ponti et 
al. 2010). On the other hand, we have found from our recent 
work that the mechanical power delivered by supernovae - 
occurring at a rate consistent with that pointed to by the 
region's current star-formation as traced by FIR emission 
— is completely sufficient to sustain the currently-observed 
non-thermal emission (~GHz radio continuum and ~TeV 
7-ray) from the ~200 pc scales of interest here. Thus, from 
the point of view proffered by the non-thermal data, it is not 
necessary that the SMBH have any significant role beyond 
the inner few pes; our investigations below confirm this in 
general. 

Finally, one of the most interesting recent discoveries in 
high energy astrophysics is of the 'Fermi Bubbles', so-called 
because these structures were revealed (Su et al. 2010) in 
~ GeV 7-ray data collected by the Fermi-LAT. The Bub- 
bles are north-south symmetric about the Galactic plane 
and centred on the Galactic nucleus. Given this morphology 
they are compellingly associated with some sort of activity 
in the GC. Given, then, the Bubbles' large angular extent 
(they rise to ±50° in b) they are enormous structures ex- 
tending 10 kpc from the plane. The Bubbles' 7-ray emission 
might be due to inverse Compton emission from a rather 
young (given short energy loss time) population of cosmic 



ray electrons. Alternatively, the emission might be due to 
hadronic collisions experienced by a hard-spectrum cosmic 
ray proton and heavier ion population (Crocker & Aharonian 
2011; Crocker 2011; Zubovas et al. 2011). We set out the 
evidence connecting the Fermi Bubbles with multi-Gyr-scale 
GC-star-formation-driven injection of cosmic ray protons into 
the Galactic halo below. 



1.1 Motivating Questions 

The preceding tour of GC and inner Galaxy phenomenology 
helps motivate a number of questions which we seek to 
address in this paper: 

(i) Gas accretion: What is the rate at which the GC 
region typically accretes gas through the Galactic plane? 
How do we understand the presence (Lubowich et al. 2000; 
Riquelme et al. 2010) of relatively pristine gas (i.e., that has 
undergone relatively little nuclear processing) in the GC? 

(ii) Star formation: What is the efficiency with which 
the GC converts gas into stars? Given the unusual conditions 
in the GC environment, is GC star-formation biased towards 
the production of massive stars? 

(iii) GC ISM conditions: Is the very hot plasma pu- 
tatively revealed by X-ray observations real or not? What 
contribution do the non-thermal ISM phases, in particular 
the cosmic rays, make to the overall energy density in the 
region? What is the dynamical importance, if any, of the 
cosmic rays? 

(iv) Outflows: There is multi-wavelength evidence (re- 
viewed below) for outfiows from the GC over size scales from 
pes to 10 kpcs. Are these outfiows all different aspects of 
the same overarching phenomenon and how are they driven? 
What is the wider importance of the GC outfiow(s) to the 
Galactic ecology? Is the material expelled from the nucleus 
lost to extra-galactic space or does it fountain back on to the 
Galactic disk? How do the recently-discovered Fermi Bubbles 
(Su et al. 2010) relate to activity in the GC? 

More broadly, we aim in this paper to produce a first draft 
of a coherent explanation of all the disparate phenomena 
listed above that is itself physically plausible and motivated. 
Overall, we shall see that it is star-formation (driven by 
secular accretion processes over long timescales) - rather 
than processes associated directly with the SMBH - that 
seems to control the overall dynamics of the GC. 



1.2 Conventions and Assumptions 

We assume a distance to the Galactic centre of 8 kpc in 
this work. We use mathematica notation: f[x] denotes / 
in a function of parameter x. Formally, the region we are 
investigating and attempting to model is that centred on 
{I, b) = (0, 0) and extending to ±0.8° in Galactic longitude 
and ±0.3° in Galactic latitude; this is the region for which 
the HESS telescope reported (Aharonian et al. 2006) a diffuse 
fiux of ~TeV 7-rays. 
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2 BACKGROUND 

2.1 Star Formation in the GC 

2.1.1 Star formation rate in the GC 

Using standard prescriptions and assuming a 'normal' initial 
mass function (IMF), ip = dN/dM, the current SFR over the 
central few degrees of the Galaxy can be estimated from the 
region's inferred Lyman continuum photon output of ~ 10^^ 
photons s~^ (Cox & Laureijs 1989) to be 0.3 — 0.6 Mo/year 
(Glisten 1989). Yusef-Zadeh et al. (2009) have determined a 
SFR over the central 400 pc that has ranged between 0.14 
and ~ 0.007 Mo/yr over the last ~ 10^ years on the basis of 
mid-infrared and other data^ These authors also estimate an 
upper limit on the SFR in the region over the last 10 Gyr of 
O.15M0/yr and claim a probable average value in the range 
0.04-0.08 Mq/jt. On the basis of 5-38 /im observations with 
the Spitzer Infrared Spectrograph of young stellar objects in 
the CMZ (and assuming a Kroupa (2001) IMF) Immer et al. 
(2011) have also recently determined a SFR in the region of 
~ 0.08 M0 /year over the last Myr. Figer et al. (2004) obtain 
a lower limit of 0.01 M©/yr on the recent GC SFR by the 
simple expedient of dividing the mass in known, recently- 
formed stars by the duration of the star formation episodes 
that formed those stars. Accounting for the fact that recent 
studies of Paschen-a emission (Mauerhan et al. 2010) have 
shown that approximately half the region's massive stars are 
located outside known clusters, a SFR of ~ O.O2M0/yr (for 
the central few tens of pes) is suggested. This SFR is close 
to the value suggested by the stellar luminosity function 
analysis of Figer et al. (2004). The detailed modelling of 
Kim et al. (2011) suggests that the SFR within the X2 orbits 
ranges between 0.04 and 0.09 Mo/year. Somewhat in tension 
with the above determinations, Yasuda et al. (2008) have 
claimed a low level of current GC star-formation activity on 
the basis of a low ratio of the [Cll] fine-structure emission line 
(due to photo-dissociation and Hll regions) to the total FIR 
emission, leading them to the conclusion that the region's 
radiative output is actually dominated by old stars. 

2.1.2 GC IMF 

Note that though the different SFR determinations listed 
above apply to regions of somewhat different sizes, given the 
GC's stellar population is highly centrally peaked (Serabyn 
& Morris 1996; Launhardt et al. 2002), there is a definite 
hint that these measures are discrepant: the region's inferred 
UV radiation output seems too high. Given that massive 
stars tend to completely dominate the production of such ra- 
diation, this discrepancy may, in fact, be an indication (Figer 
et al. 2004) that the region's star formation is biased towards 
the production of massive stars (Morris 1993; Maness et al. 
2007), consistent with independent indicators (Figer et al. 

1 The 0.007 Mq/jt reported by Yusef-Zadeh et al. (2009) is 
obtained from their estimate of the 1.4 GHz radio continuum flux; 
such an estimate suffers from the problem that cosmic ray electrons 
are removed from this system before they can radiate much of 
their energy: see below and also Crocker et al. (2011a, b). Put 
another way, were this advective removal of cosmic rays electrons 
not taking place one would derive a RC estimate of the SFR much 
closer to determinations relying on other data. 



1999) that the region's initial mass function is significantly 
flattened. That star-formation in the GC be biased towards 
the formation of more massive stars is a rather natural pre- 
diction: given the environment, any or all of the region's 
strong tidal forces, high gas pressures, and magnetic fields 
might be expected (Morris 1993; Lis et al. 2001) to signifi- 
cantly alter the dynamics of the collapse of molecular gas into 
stars. Whether this expectation is supported observationally 
remains, however, a topic of hot debate. Indeed, the GC, as 
the site of some of the most active massive star-formation in 
the Galaxy, has been a natural battle ground in the debate 
over whether the IMF is truly universal (e.g., Bastian et al. 
2010; Lockmann et al. 2010) or is flattened (or has a higher 
low mass cut-off to the formed stellar population) in the GC 
(Morris 1993; Figer et al. 1999; Figer et al. 2004; Maness et 
al. 2007) and other star-burst-like environments. 



2.1.3 Continuous and steady star formation in the GC 

Importantly for our purposes, the luminosity function analy- 
sis of Figer et al. (2004) also favours a SFR that has been 
sustained at more-or-less the current value for a timescale 
approaching 10 Gyr. In fact, there seems to be accumulat- 
ing evidence from different directions (Maness et al. 2007; 
Lockmaim et al. 2010; Kim et al. 2011) that the GC has been 
continuously forming stars over this sort of timeframe. Indeed, 
we (Crocker & Aharonian 2011) have recently suggested on 
the basis of our modelling of the non-thermal emission Fermi 
Bubbles (see below) that the currently-observed star forma- 
tion rate in the GC is typical of the system's time-averaged 
value over the last ~ 8 Gyr. 

A corollary of this sort of picture is that the drama 
associated with most of the accumulation of the mass of 
the SMBH at Sgr A* is pushed back to highish redshifts, 
probably accompanied by the formation of most of the long- 
lived stellar population of the bulge. In this context, the fact 
(Gilmore et al. 2002) that the Milky Way underwent its last 
major merger activity ^11 Gyr ago and has subsequently 
experienced rather quiescent (Yin et al. 2009) evolution is 
significant. On the other hand, recent modelling (Purcell 
et al. 2011) shows that at least some of the continuing gas 
feeding to the inner Galaxy could ultimately be driven by 
the on-going, minor merger activity the Galaxy experiences 

Other independent evidence that the GC SFR has been 
steady over long timescales comes from the observation (dis- 
cussed further below: §2.3.6) of separate parcels of highly- 
ionized high-velocity gas (in UV absorption spectroscopy 
along the sight lines to distant quasars) at low Galactic lon- 
gitude but varying Galactic height (both north and south 
of the plane), both emerging from and, apparently, falling 
back on to the GC (Keeney et al. 2006). Significantly, these 
gas parcels apparently form part of a Galactic fountain and 
can be inferred to have reached (or will reach) the same 
maximum height from the plane of 12 ± 1 kpc but must have 
been launched at different times in the past ranging from 
20-50 Myr to more than 800 Myr, suggesting the operation 
of a common launching mechanism over at least the latter 
timeframe. The most recent UV measurements (towards a 
post-AGB in the inner Galaxy; Zech et al. 2008) suggest, 
moreover, that at least some of this material is of super-solar 
metallicity. This suggest both an inner-Galaxy origin and 
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may point to a star-formation - rather than AGN - origin 
to the outflow(s) (cf. Su ct al. 2010). 

2.2 X-ray observations 

In its plasma phase, X-ray continuum and Fe line observations 
apparently reveal a two-temperature plasma containing 'hot' 
(~ 1 keV) and (mysteriously) 'very hot' (6-9 keV) components 
(Koyama et al. 1989; Yamauchi et al. 1990; Kaneda et al. 
1997; Muno et al. 2004; Belanger et al. 2004). The X-ray 
emission from the putative very hot component is strongly 
concentrated within the inner ~150 pc (in diameter) of the 
Galaxy (Yamauchi et al. 1990; Belmont et al. 2005). As first 
observed by Spergel & Blitz (1992), there may be pressure 
equihbrium (at (3 - 6) x 10^ K cm"^: Koyama et al. 1996; 
Muno et al. 2004) between the kinetic pressure of the putative 
very hot plasma phase and the virial pressure implied by the 
turbulent motions of the molecular gas. 

Prima facie, the very-hot plasma presents a severe ener- 
getics problem, however: (assuming it is a hydrogen plasma) 
its sound speed at ~1500 km/s would be considerably in 
excess of the local escape velocity of ~ 900 km/s (Muno et 
al. 2004) suggesting it should escape (Yamauchi et al. 1990) 
on a short timescale. This suggests a steady-state situation 
would require a power considerably in excess of 10"*° erg/s 
to sustain the outflow. A second difficulty is that there is no 
widely-accepted mechanism to heat the plasma to more than 
a few keV; Galactic disk supernova remnants, in particular, 
do not seem to heat plasma beyond ~ 3 keV after a couple of 
centuries (hotter temperatures at earlier times are possible 
but the smooth distribution and overall energy of the GC's 
putative hot plasma cannot be reconciled with such a young 
explosion: Muno et al. 2004; Belmont & Tagger 2006). 

There is no universally-accepted resolution to these 
anomalies. One interesting suggestion is that the 8 keV emis- 
sion is due to a very hot helium plasma which would be 
gravitationally bound (Belmont et al. 2005). Another sug- 
gestion is that the 'plasma' is illusory, the emission actually 
being attributable to unresolved point sources(Wang et al. 
2002). Recent, deep Chandra observations around I = 0.08, 
b =1.42 (taken to be typical of the so-called X-ray Ridge) 
support this sort of picture (Revnivtsev et al. 2009). On the 
other hand, the situation within the inner ~150 pc - where 
the 6.7 keV Fe line emission strongly peaks (Yamauchi et al. 
1990) - may be quite different to that pertaining elsewhere in 
the Galaxy (Dogiel et al. 2009a). A deep observation (Muno 
et al. 2004) of the inner 17' with Chandra could only ex- 
plain ^ 40% of the X-ray flux as due to dim point sources. 
Moreover, recent results obtained with the SUZAKU X-ray 
telescope continue to clearly suggest (Koyama et al. 2007; 
Dogiel et al. 2010; Koyama 2011) the existence of a hot 
plasma covering at least the central 20'; this issue, therefore, 
has remained unresolved. Below the modelling we present 
shows how the GC's star-formation activity might be able to 
sustain such a plasma in steady-state. 



2.3 Evidence for a Galactic Centre Outflow 

Observationally, there is evidence on multiple scales and at 
many different wavelengths for an outflow or outflows from 
the GC, some of which we review briefly below. The idea 



of an outflow has also received theoretical support from our 
recent work (Crocker et al. 2011a; Crocker & Aharonian 
2011; Crocker et al. 2011b) which we also briefly review. 

2.3.1 GC Lobe and CMZ molecular halo 

Radio and optical recombination line observation observa- 
tions by Law et al. (2009) reveal 2 x 10^ Mg of warm (~ 4000 
K), ionized gas extending up to ~ 1° north from the plane 
above the GC. This gas is nested within a shell formed by the 
so-called Galactic Centre Lobe (GCL; Sofue & Handa 1984) 
detected as a non-thermal radio continuum source at ~GHz 
frequencies (Law 2010, and references therein) and visible 
up to at least 10.5 GHz (Sofue 1996; Crocker et al. 2010a; 
refer to the radio continuum contours in fig. 1 from Pohl 
et al. 1992: the GCL rises between Z ~ 0°.2 and I ~ -0°.7, 
i.e., roughly above the Radio Arc and Sagittarius C). An 
outermost shell around this structure of entrained dust and 
PAHs is detected at MIR wavelengths (Bland-Hawthorn & 
Cohen 2003), some emission revealing helical topologies pre- 
sumably tracing a complex magnetic field structure (Morris 
et al. 2006). Finally, CO and CS line emission from the region 
reveals molecular gas extended along spurs north (Uchida 
et al. 1994) and south of the plane and more-or-less coinci- 
dent with (actually slightly inside-of) the radio continuum 
features and apparently rotating (Sofue 1996). The mass of 
such molecular gas in the GCL has been estimated to be at 
least ~ 3 X 10^ A'/©. 

The number density of the warm gas from the radio 
recombination line observations can be estimated to be ~ 10^ 
cm"'' and its pressure P/ks ~ 7 x 10® K cm"'' would put it 
in or close to pressure equilibrium with the other GC ISM 
phases (including the very hot plasma were it real; Law et 
al. 2009). The gas is also of high metallicity and parts of 
the lobe's RC emission suffer from Hi absorption (Law et 
al. 2009); both these factors clearly point to the structure's 
location in the GC. The warm gas has, however, a low filling 
factor, / ~ 10~*, much smaller than typical for this phase 
in the Galactic disk but interestingly comparable to that 
inferred for other Galactic outfiows (Law et al. 2009, and 
references therein). 

Radio continuum observations also lend support to the 
notion that the GCL represents an outflow (Law 2010): the 
non-thermal spectrum of the GCL steepens as a function of 
increasing Galactic latitude, a clear sign of a synchrotron- 
emitting electron population that is ageing as it is transported 
from the plane (cf. the recent work by Heesen et al. 2009, 
on the star-burst system NGC 253). 

We can use the estimated ~ 2 x 10^ Mq of warm, ionized 
gas fllling the GCL (Law 2010) to a height h ~140 pc north 
of the GC to arrive at a lower limit to the mass flux, M in 
an outflow: with Mqcl ^ Mh/v^M > Mobs = 2 x 10^ Mq, 
we find i\f ^ 0.3 A'/g/year xvmind/WO km/s (assuming a 
similar, but so-far unobserved, distribution of ionized gas 
south of the plane from the GC - though see §2.3.2)^^. 

The above is likely a conservative estimate as it neglects 

^ Scaling the results of Martin (2005, figure 6) according to the 
GC's estimated SFR areal density, we find that the expectation 
afforded by observations of external galaxies is that the GC should 
drive an outflow with a speed of '~400 km/s. 
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Figure 1. Colours: ~ TeV map of the GC region (Aharonian et 
al. 2006) after removal of point-like sources coincident with Sgr A* 
and the SNR GO. 9+0.1 and (white contours) CS(l->0) molecular 
line emission (a tracer of the total H2 coliunn); black contours: 2.7 
GHz unsharp-masked radio continuum image with 9.4' resolution 
(Pohl et al. 1992). The radio continuum spurs above I ~ 0°.l (the 
Radio Arc) and ~ — 0°.7 (Sgr C) define the non-thermal GC Lobe 
structure. Note that the footprints of the radio continuum spurs 
are coincident with regions of high molecular column and ^TeV 
7-ray emission and star-formation. HESS data figure is adapted 
by permission from Macmillan Publishers Ltd: [Nature, 439, 695], 
copyright (2006). For 2.7 GHz radio contours credit: Pohl et al. 
(1992), A&A, 262, 441, 1992, reproduced with permission ©ESQ. 

mass in other gas phases. In fact recent CO(2 — >■ 1) and 
CO(l — >■ 0) line observations with the Nanten-II telescope 
point to a halo of molecular gas around the entire CMZ with 
height ~ 1° (or 100-200 pc) and total mass few x lO'^ Mq 
(Yasuo Fukui, private communication). Spitzer IRS spectra 
also reveal high-latitude 17 and 28 /xm H2 emission lines 
around the CMZ, consistent with such a molecular gas halo 
(Mark Morris, private communication). Finally, older OH 
absorption line observations of the region (Boyce & Cohen 
1994) also point to the existence of a ~ 200 pc molecular 
halo around the GC and reveal individual molecular spur 
features coincident with radio continuum features (including 
the GCL). 

2.3.2 GC 'worms' 

Koo et al. (1992) made an inventory of Galactic 'worms': 



isolated structures appearing in both Hi and far infrared (60 
and 100 /xm) and generally identified as the walls around 
superbubbles that have broken through the thin gaseous disk 
of the Galaxy (to both north and south). Many of these 
structures are coincident with Hll regions and 408 MHz RC 
emission. Koo et al. (1992) list four worms coincident with 
the GC region under consideration in this paper, with three 
(GW 357.4 + 5.4, 358.0 - 5.3, 0.5 - 5.9) appearing around 
vlsr = km/s ("LSR" denotes local standard of rest). The 
inferred mass of GW 0.5 - 5.9 (which projects south of the 
Galactic plane) is 2 x 10^ Mq , interestingly close to the mass 
of warm ionised gas observed in the GC lobe (which projects 
north) . 

2.3.3 Extended molecular gas and dust distribution 

There is multi-wavelength evidence for extended - and proba- 
bly outflowing - cold gas around the GC. Surveys of molecular 
lines and dust continuum emission at mm, sub-mm, and mid- 
infrared wavelengths (Oka et al. 1998; Tsuboi & Miyazaki 
1998; Pierce-Price et al. 2000; Bland-Hawthorn & Cohen 
2003; Stolovy et al. 2006; Oka et al. 2010) reveal filaments, 
arcs, and shells, indicative of local, turbulent sources (Tanaka 
et al. 2007) and explosive events. A large-scale CO J=3-2 
survey of the CMZ with the Atacama Submillimeter-wave 
Telescope Experiment (Oka et al. 2010) detects an unusual 
population of high- velocity {Av ^ 50 km/s) compact clouds 
(HVCCs) in the GC. 24 micron observations have revealed 
dust emission from a fascinating 'double helix' structure 
(Morris et al. 2006) ~ 0.°7 degrees north of the GC which 
has recently been shown^ to exhibit kinematically-related 
molecular emission and linearly-polarized radio continuum 
emission consistent with a highly-ordered magnetic field 
structure(Tsuboi & Handa 2010); ^^CO (1 0) molecular 
line emission observations demonstrate the existence of giant 
molecular loops (GMLs) with large velocity dispersions in the 
region, argued to originate in the phenomenon of magnetic 
fiotation controlled by the Parker (1966) instability (Fukui 
et al. 2006). Importantly for our purposes here, the most 
recent ^^CO (1 — >■ 0) and (2 — >■ 1) molecular line observations 
by the Nanten-II group reveal an extended halo of molecular 
gas around the entire CMZ with total mass of few x 10® Mq 
of H2, with more in individually-identifiable outflows'*. 

2.3.4 Evidence for disk-halo connection m the GC 

Another interesting piece of the puzzle of GC gas dynamics 
emerges from studies indicating the presence of relatively pris- 
tine (i.e., relatively nuclear-unprocessed) or even primaeval 
gas in the GC. A first piece of evidence suggesting accretion 
of such gas into the region comes from observation of the 
J = 1^0and2->1 lines of DCN in the '50 km s"^' 
molecular cloud located ~10 pc from the true GC (Lubowich 
et al. 2000). Deuterium is destroyed in stellar interiors and 

Rei Enokiya, talk presented at 'The emerging, multi-wavelength 
view of the Galactic Centre Environment', Heidelberg, Germany, 
October 17-20, 2011. 

^ Yasuo Fukui, private communication, and Kazufumi Torii, talk 
presented at 'The emerging, multi-wavelength view of the Galactic 
Centre Environment', Heidelberg, Germany, October 17-20, 2011. 
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- unless it is continuously replenished in a heavily astrated 
region like the GC - its abundance relative to H should be 
at the few parts per trillion level; Lubowich et al. (2000) 
determined levels five orders of magnitude higher than this. 

This inference of the infall of fresh gas into the GC has 
recently found interesting confirmation in determinations of 
high ^^C/^^C isotopic ratios (through measurement of the 
1 lines of HCO+, HCN, HNC and their "C isotopo- 
logues) along a number of GC sightlines. Overall the ^^C/^''C 
ratio is known to exhibit a gradient from high (80-90) to 
low (20-25) values going from the galactocentric distance of 
the solar system to the GC (Wilson 1999). This is consis- 
tent with the general picture that gas in the outer Galaxy 
should have experienced less nuclear processing than that 
in the inner: ^^C is formed in first-generation, metal-poor 
stars over shortish timescales whereas ^''C is formed via CNO 
processing by low or intermediate mass stars of ^^C seeded 
by earlier stellar generations (Riquelme et al. 2010, and refer- 
ences therein). Riquelme et al. (2010) have determined lower 
limits to ^^C/^'^C isotopic ratios for a number of places in the 
CMZ that are close to the values measured in the local ISM 
and inconsistent with the general gradient identified above. 
The implication of this - as for the deuterium observations - 
is that some relatively less-processed gas is, somehow, finding 
its way directly to the GC. 

Riquelme et al. (2010) have found evidence for such 
unprocessed gas, in particular, in the footprints of a number 
of GMLs identified by the Fukui group and, in general, in 
gas whose phase space position places it in the XI orbit 
family or in transit from XI orbits to X2 (recent higher 
resolution data from Jones et al. 2011, seems to confirm 
this latter finding; see their fig. 18). Indeed, the measured 
i2(-i^13q jgotopic ratios may be more consistent with the idea 
(Morris 2006; Torii et al. 2010; Riquelme et al. 2010) that the 
shocks induced in the rising portions of the loops sweep up 
and compress rarefied atomic gas in the halo above the GC, 
leading to rapid cooling and condensation into molecular gas. 
Whatever the mechanism, the presence of relatively pristine 
gas in the GC, represents compelling evidence that some gas 
is being rather directly accreted out of the halo on to the 
GC. 



2.3.5 Hi 'outflow(s)' from GC 

Making sense of 21 cm line data towards the GC is difficult. 
The interpretation of features in such data as indicating some 
sort of outfiow or multiple ejections from the nucleus has a 
venerable history (e.g., Sanders & Wrixon 1972; Mirabel & 
Franco 1976). Much of the apparently anomalous kinematics 
of the individual atomic hydrogen features within jbj ^2° 
over the |Z| ^ 1° longitude range of relevance here has subse- 
quently been reinterpreted as motion governed by the inner 
Galaxy's tilted Hi disk identified by Burton & Liszt (1978). 
Even this latter work, however, found tentative evidence for 
streams of material moving out along the polar axis of this 
disk (tilted by 12° with respect to the vertical in the model) 
away from the nucleus at ~ 200 km/s and representing a 
total mass ^ 10^ Mq. We remark in passing that these op- 
posing streams - moving into the north east and south west 



quadrants'' - are roughly aligned to the corresponding, inner 
edges of the Fermi Bubbles (see below, §2.3.9) and the bicon- 
ical X-ray feature identified in ROSAT data by Sofue (2000) 
(see below, §2.3.8); Sanders & Wrixon (1972) identified an 
Hi feature on similar size scales that seems to be similarly 
coincident with the noith-west edge of the northern Fermi 
Bubble and the same X-ray structure. 

On much larger scales, individual, high and intermedi- 
ate velocity Hi clouds or cloud complexes are seen all over 
the sky (e.g., Winkel et al. 2011, and references therein); 
the distance to these features - and their inferred masses 
and sizes - is notoriously difficult to pin down. A massive 
literature describing these data exists; it is sufficient for 
current purposes to mention that a small fraction of such 
clouds may be associated with material ejected by a GC 
fountain. The recent analysis of Winkel et al. (2011) of the 
high-velocity cloud complex Galactic center negative (cov- 
ering the south-east complex out ~ 80° from the GC) finds 
a distinct sub-population of clouds which, while separated 
from the GC by up to 80° , seems to be kinematically related 
to high-velocity dispersion gas in the central few degrees of 
the Galaxy; this might be tentatively ascribed to a large- 
scale GC ejection event or outfiow, though analysis of this 
structure is on-going. 

2.3.6 UV absorption spectroscopy: evidence for a GC 
fountain 

As prefaced above, UV spectroscopic data suggests both ejec- 
tion of highly-ionized material from the GC (Bland-Hawthorn 
& Cohen 2003) and the subsequent fountaining-back of some 
of this material (Keeney et al. 2006). At least some of this 
material seems to have super-solar metaliicity (Zech et al. 
2008). Consistent with this evidence, we will see below that 
our modelling naturally predicts an outfiow speed less than 
the gravitational escape speed. 

2.3.7 GC spur 

On the basis of an analysis of unsharp-masked 408 and 1408 
MHz radio continuum data, Sofue et al. (1989) identified 
a large, non-thermal (but hard spectrum) radio feature ex- 
tending north of the GC up to latitudes of ~ 20° that they 
labelled the Galactic Centre Spur (GCS). Jones et al. (2012) 
have recently shown that this same feature is visible in po- 
larised emission in WMAP microwave data (at 23 and 33 
GHz). There is no obvious, corresponding feature seen to 
Galactic south in radio or microwave frequencies. Rather 
interestingly in light of the discovery discussed in the next 
sub-section, by eye the GCS appears also partially coincident 
with a 7-ray feature extending north from the GC that was 
claimed previously on the basis of EGRET data (Hartmann 
et al. 1997, though note the statistical significance of this 
feature could not be established). 

Also notable is that - despite its length and curvature - 
the GCS remains well-coUimated over its length having an 
almost constant width of 1.5 — 2° (Sofue et al. 1989). This 
essentially matches the width of the GC star-forming region 

The south-west feature is tracked by Mirabel & Franco (1976) 
out to 6 = -22°. 
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we model in this paper. Exactly how the structure remains 
coUimated over distances of ^ 4 kpc is mysterious; regardless 
we make the point that it apparently represents a channel 
for the delivery of non-thermal particles to large distances 
from the plane with minimal adiabatic energy losses. 

At 1.4 GHz the GCS divides into two strong radio spurs 
which very plausibly - but not definitely - join on to the 
radio spurs seen (Pohl et al. 1992) above Sgr B and Sgr C at 
2.7 GHz (see fig. 1). In addition to the above, a number of 
other pieces of evidence mark the GCS as a unique feature 
and suggests its GC location: 

(i) the feature terminates in the Galactic plane; 

(ii) it is the brightest radio continuum spur after the North 
Polar Spur (Sofue et al. 1989), likely a local ISM feature; 

(iii) while polarised emission from this structure is de- 
tectable at microwave frequencies (Jones et al. 2012), po- 
larised emission from the GCS disappears by 1.4 GHz: such 
behaviour is consistent with the 'magnetic horizon' effect, i.e., 
the Faraday depolarisation due to the turbulent ISM which 
renders polarised emission at ~GHz undetectable beyond a 
few kpc through the GC plane (thus the GCS is likely at at 
least this distance). 

Finally, the GCS exhibits noticeable curvature to Galac- 
tic west; this curvature can be coherently explained within 
the general idea that the feature is due to a rather slowly 
moving, star- formation-driven outflow (~ 1000 km/s). In 
particular, if the feature is associated with an individual 
star-formation 'event' occurring in the central ~100 pc star- 
forming gas ring, rotating at ~ 100 km/s, (Molinari et al. 
2011), then differential rotation would indicate a formation 
timescale ~ tt 100 pc/100 km/s ~ 3 x 10^ year (note that 
the inferred outflowing corkscrew has not yet executed one 
full turn) and an outflow speed ~ 900 km/s (cf. Heesen et al. 
2011, on the helical magnetic field structure of the nuclear 
outflow from NGC 253). These inferred parameters seem 
eminently reasonable; further study might be able associate 
the feature to the formation of a particular, GC super stellar 
cluster. Note that while the inferred GCS outflow speed is 
somewhat faster than for the general outflow we identify (see 
below), the speed is rather well matched to that required 
to generate the sort of high latitude, high metallicity, and 
highly-ionized gas features identified by UV spectroscopy 
(Keeney et al. 2006; Zech et al. 2008). 

2.3.8 X-ray evidence of a giant, GC-centred bicomcal 
structure 

Sofue (2000) identified a biconical structure extending north 
and south from the GC in 1.5 keV ROSAT data that extends 
out to |6| ~ 20*^ above and below the plane. The dynam- 
ics of this structure were extensively investigated by both 
Sofue (2000) and later Bland-Hawthorn & Cohen (2003). So- 
fue (2000) suggested an association with very-large-angular 
scale structures in the 408 MHz sky, in particular, the North 
Polar Spur and suggested these multi-wavelength features 
were related to a nuclear star-burst that released ~ 10^^ 
erg energy, filling out a giant hyper shell centred on the GC. 
Bland-Hawthorn & Cohen (2003) similarly drew a connection 
between the biconical X-ray feature and multi-wavelength 
data on different scales including the MIR features mentioned 
above and again suggested a likely star-burst origin for the 



apparent outflow, though with somewhat more modest ener- 
getics (~ 10^^ erg). 



2.3.9 Fermi Bubbles 

One of the most interesting recent discoveries in high-energy 
astronomy is of the 'Fermi Bubbles' (Su et al. 2010; Dobler et 
al. 2010) introduced earlier. The Bubbles are characterised by 
a rather uniform intensity and an unusually-hard spectrum, 
dF^i/dE-y oc E:^'^'^ and have a total luminosity 4 x 10'^^ erg/s. 

Many researchers have focused on the general idea that 
the 7-ray emission from the Fermi Bubbles arises from the 
inverse-Compton (IC) emission from a (mysterious) popula- 
tion of cosmic ray electrons. Given that the spectrum of the 
Bubbles displays no obvious variation with Galactic latitude, 
however, it is necessary that the photon background being 
up-scattered by this putative electron population is the CMB. 
This, in turn, implies that the electrons have an energy scale 
~TeV and consequently short IC loss times, ~ 10® years 
(Crocker & Aharonian 2011). Given the vast extension of 
the Bubbles (~ 10 kpc from the plane), these electrons ei- 
ther have to be delivered very quickly - presumably on an 
AGN-type outflow originating at Sgr A* (Guo & Mathews 
2011) - or accelerated in-situ by first (Cheng et al. 2011) or 
second-order (Mertsch & Sarkar 2011) Fermi acceleration 
processes. 

We (Crocker & Aharonian 2011) have recently consid- 
ered the alternative explanation that the Bubbles' 7-ray 
emission arises from the hadronic collisions of a population 
of cosmic ray protons (and heavier ions) populating their inte- 
riors. Because of the long proton loss times on the low-density 
plasma of the Bubbles this escapes the timing difficulties 
facing any leptonic mechanism. Our explanation requires i) 
(given adiabatic and ionisation energy losses) a total cosmic 
ray hadron power ~ 10"^^ erg/s that ii) (essentially because 
of the same long loss time referred to above and the conse- 
quently long time required to reach steady state) has been 
injected quasi-continuously into the Bubbles for a timescale 
of 8 Gyr. These requirements are precisely matched by the 
GC CR outflow that we (Crocker et al. 2011b) identified on 
the basis of completely independent considerations to do with 
observations at radio continuum and TeV 7-ray wavebands 
of the inner ~ 200 pc of the Galaxy. This putative solution 
fits nicely from a number of other perspectives: 

(i) The hard-spectrum of the emission is also explained: 
by construction, the cosmic rays injected into the Bubbles 
are trapped so there is no energy-dependent escape process 
acting to modify the in-situ, steady state distribution away 
from the injection spectrum and the daughter 7-rays will 
trace this hard, parent proton distribution. 

(ii) On the other hand, vr^-decay kinematics enforces a 
down-turn below ~GeV on a spectral energy distribution 
plot of the emitted 7-radiation; such a downturn is robustly 
detected, at least qualitatively, in the Bubbles' spectra (Su 
et al. 2010). 

(iii) The total enthalpy - the energy required to supply 
the final internal energy of the Bubble and do the p dV 
work of infiation against the pressure of the external medium 
- of a slowly infiated Bubble is (e.g., Hinton et al. 2007) 
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given hy H = 7/(7 - 1) p V ~ (2 - 4) x 10^^ erg^ 
(where 7 = 4/3 — 5/3 and we adopt 4 x 10'* K cm""^ as 
the pressure in the Galactic disk towards the inner Galaxy: 
Kasparova & Zasov 2008, and the total volume of the Bubbles 
is ~ 2 X 10^^ cm~^). Setting Eout tinf = H where Eout is 
the rate at which the outflow does P dV work and tinf is 
the inflation timescale, we flnd ti^f ~10 Gyr (i5out/10^^ 
erg/s )~^\ where we normalise to the typical power we flnd 
that the GC delivers into freshly-accelerated cosmic rays 
(see below). In reality, the cosmic rays, advected, frozen- 
in magnetic fleld, and injected plasma will all contribute 
to inflating the Bubbles but the important point is that 
this inflation will happen over the same multi-Gyr timescales 
independently suggested hy considerations around reproducing 
the j-ray phenomenology of the Bubbles given the long pp 
timescale. (While the outflow is likely to have a total power 
> 10*° erg/s, both cosmic rays and the plasma injected 
into the outflow will experience significant radiative losses 
over these long timescales so adiabaticity is not satisfied: 
see below.) Note that the alternative postulate that the 
Bubbles are inflated by much shorter duration and much 
higher mechanical luminosity output from a recent AGN 
phase of the SMBH comes close to violating a limit on such 
activity claimed (Lubowich et al. 2000) on the basis of the 
DCN line observations referred to above. 

(iv) The total plasma mass of the Bubbles is 10* Mq 
(Su et al. 2010) - this mass can also be explained given the 
rate of mass flux in the GC outflow and assuming the same 
long timescales (as we conflrm in detail below). Also note 
that the total power fed into the base of the Bubbles by the 
outflow can sustain the thermal (X-ray) radiation from the 
Bubbles (Crocker & Aharonian 2011). 

(v) Dynamically, the Bubbles end up being slightly over- 
pressured but slightly under-dense with respect to the sur- 
rounding halo plasma, with internal energy density supplied 
approximately equally by cosmic rays and their interior hot 
plasma. They can, therefore, be expected to rise slowly under 
buoyancy. 

(vi) The hadronic scenario naturally predicts concomitant 
secondary electron production within the Bubbles; these sec- 
ondaries would synchrotron-radiate on the Bubble's magnetic 
fleld, thereby explaining the coincident (at lower Galactic 
latitude) 'WMAP haze' detected (Finkbeiner et al. 2004; 
Dobler & Finkbeiner 2008) at microwave frequencies. 

Of course, all this requires that the Fermi Bubbles are 
very old structures - almost as old as the Galaxy - and that 
they can trap TeV cosmic rays for multi-Gyr timescales. In 
fact, our scenario implies that they would be calorimeters 
for GC activity over the history of the Milky Way. This is 
an interesting prospect indeed. 

We finally remark on a very recent development: 
Finkbeiner and co-workers^ have recently claimed the discov- 
ery of 7-ray substructure within the Bubbles, in particular 
two counter-propagating jet-like features intersecting the GC 



^ Note that Guo & Mathews 2011 find from their modelling that 
typical total AGN energy release of ~ 10^^ erg is required in a 
fast-inflation scenario with AGN jets. 

^ Douglas Finkbeiner, talk delivered at The Emerging, Multi- 
Wavelength View of the Galactic Centre Environment, Heidelberg, 
Germany, October 2011. 



and extending into both Galactic hemispheres (slightly west 
of north and east of south, respectively) with even harder 
spectra than the overall spectrum of the Bubbles. By eye, 
these features seem at least partially coincident with the 
GCS feature identified in radio continuum and polarised 
microwave emission and discussed above; in contrast to the 
widescale Bubble emission, they may originate in IC emission 
from young, primary electrons carried out of the GC system. 

2.3.10 Non-thermal evidence for an outflow 

As briefly reviewed above, the GC displays extended, diffuse 
TeV emission (Aharonian et al. 2006) spatially correlated 
with the column of molecular gas over the central ~ 1.5° in 
Galactic longitude. On even wider scales than for the TeV 
emission (~ 6° in Z, ~ 2° in b: LaRosa et al. 2005; Crocker et 
al. 2010a), radio continuum observations show that the GC 
is a distinct source of diffuse, ~GHz, non-thermal emission 
(of which the GC lobe forms a part). Such emission must be 
due to the synchrotron losses experienced by a wide-spread 
population of cosmic ray electrons inhabiting the GC. 

Despite the fact of this wide-spread non-thermal emis- 
sion, the GC is actually significantly underluminous in both 
radio continuum and ~TeV (and ~GeV) 7-ray wavebands 
given the amount of star-formation currently going on there - 
as we now explain. Firstly, placing the GC on a plot of its 60 
^m vs. 1.4 GHz luminosity, one determines that radio contin- 
uum emission from this system falls one order of magnitude 
(i.e., ~ 4(t) short with respect to the expectation afforded 
by the FIR-radio continuum correlation (e.g., Condon 1992). 
Equally, confronted with the theoretical expectation for the 
numerical scaling (Thompson et al. 2007) between a star- 
forming galaxy or region's (non-thermal) 7-ray and its FIR 
emission (expected were the region calorimetric to the accel- 
erated cosmic rays), the GC is significantly in deficit. In fact, 
its TeV luminosity is at the level of ~ 1% expectation (the 
system's GeV emission, as measured by the Fermi satellite 
(Chernyakova et al. 2010), is at about 10% of expectation, 
but substantially polluted by point sources in the field). 

As we have discussed at length previously (Crocker et 
al. 2011a,b), the explanation for these non-thermal deficits 
is that the GC is not a calorimeter for either the cosmic ray 
protons or electron populations it accelerates: some sort of 
transport process is acting to remove the non-thermal parti- 
cles quickly enough that they do not have the opportunity to 
lose their energy radiatively in situ. This explanation seems 
reasonable - and, indeed, given the small scale of the GC re- 
gion under consideration, we are probing length scales where 
the FIR-RC correlation is known to break down (Hughes 
et al. 2006) - but the 'usual' explanation that cosmic ray 
diffusion is removing the particles is not available to us. This 
follows for the simple reason that - given the hardness of 
the detected non-thermal radiations - there is no evidence 
for diffusion steepening of the emitting non-thermal particle 
populations (cf. the situation in the Galactic plane). Indeed, 
the steady state particle populations appear to have a spec- 
trum - dN/dE oc or so - completely consistent with 
the expectation for that at injection (following first-order 
Fermi acceleration at astrophysical shocks) . Thus, if some 
process is acting to transport particles away - as apparently 
required on the basis of the evidence described above - this 
process must act without prejudice as to particle energy. 
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This requirement is naturally met by a large-scale outflow or 
wind of a few hundred km/s. Note also that the implication 
of this reasoning is that the GC is losing ~ 10^® erg/s in 
non-thermal particles to the outflow: this is precisely the 
power required to sustain the 7-ray emission from the Fermi 
Bubbles in steady state as identified above (and also to in- 
flate them over the same ~ 10^° year timescale required for 
the pp-in-saturation explanation of the 7-rays). 



3 MODELLING 

Given all the forgoing, we seek to understand the non-thermal 
~GHz radio continuum and ~TeV 7-ray emission detected 
from the GC region. To this end we have developed a single- 
zone model of the injection, cooling, and escape of relativistic 
particles from the region. Cooling and escape processes - con- 
trolled by the environmental conditions as described by our 
parameter space - form the steady-state, non-thermal parti- 
cle populations we model; we also make the theoretically and 
empirically-motivated assumption that cosmic ray popula- 
tions are injected into the ISM by their accelerators as power- 
laws in momentum. Relevant cooling processes are hadronic 
(pp) collisions and ionizing/Coulomb collisions for high- and 
low-energy protons, respectively, and ionizing/Coulomb col- 
lisions, bremsstrahlung, synchrotron, inverse Compton (IC) 
emission for electrons, and adiabatic losses for all particle 
types. Our code also accounts for advective particle loss. We 
track production of electrons and positrons ('secondary elec- 
trons') through charged meson decay following pp collisions 
and the radiation they produce. Relevant radiative processes 
are, at radio continuum (and microwave) wavelengths, syn- 
chrotron emission by electrons and, at 7-ray wavelengths, 
bremsstrahlung and IC, by electrons and neutral meson decay 
following pp collisions for protons. 

Self-consistently, the radiation from the modelled steady- 
state non-thermal particle populations should reproduce the 
emission we detect from the GC region*; we search over our 
parameter space with a minimisation procedure (instanti- 
ated within mathematica), trying to reproduce the particle 
populations and environmental conditions that do this. Our 
modelling of the non-thermal particles - in principle, cosmic 
ray protons and heavier ions and electrons, but here taken 
to be simply protons and electrons for simplicity - their radi- 
ation, and the secondary particles they produce, is largely as 
described in (Crocker et al. 2011a,b) with some extensions 
and modifications as described below. 

3.1 Assumptions 

One difference with our previous modelling is that - in order 
to cut-down the size of the parameter space we need to 
search - we assume that the ratio of electrons to protons at 
injection, K^p, follows the theoretical expectation (Bell 1978) 
for momentum power laws with identical spectral index 7e = 
7p = 7 and assuming equal overall numbers of electrons and 
protons accelerated into non-thermal populations (e.g., for a 
typical best-fit spectral index value 7 — 2.2 this expectation 

* See Crocker et al. (2011b) for modelled broad-band spectra of 
the region. 



IS Kep 

[7 = 2.2] ~ 0.01 at TeV). This simpfification is justified 
given that we found in our previous modelling (Crocker et al. 
2011b) that the best-fit values of a fioating Kep correspond 
exactly to this theoretical expectation. 

Other important assumptions of our modelling of the 
non-thermal particle populations are that the system is in 
quasi-steady-state; this condition was previously shown by 
us to be empirically reasonable (refer fig. 4 of Crocker et al. 
2011b). Essentially this is guaranteed if the time between 
injection events is the smallest relevant timescale, shorter, in 
particular, than particle cooling and escape timescales. This 
condition is met over the parameter space of relevance with 
possible exceptions for the case of very high energy electrons 
(> TeV) in strongish magnetic fields (> 100 pG) and where, 
in addition, rare hypernovae {Esn 10^^ erg as defined 
below) make a significant contribution to the total, time- 
averaged energy budget of the region. Practically, however, 
this case is not a concern given that we find 7-radiation (IC 
emission) from these high-energy electrons is subdominant to 
hadronic 7-ray emission over best-fit regions of the parameter 
space. Note that, even if the steady state condition is not 
formally met, the condition is too stringent as it neglects 
the finite timescale for which each hypernova explosion is 
an effective accelerator. Indeed, it is possible (Crocker et al. 
2011b) - some have even argued likely (Melia & Fatuzzo 
2011; Amano et al. 2011) - that the unusually turbulent 
ISM conditions in the ISM mean that particle acceleration 
occurs not (or not exclusively) in association with individ- 
ual SNRs per se but rather via second order (stochastic) 
acceleration on diffuse ISM turbulence. If this is the case, 
the region should better be regarded as a single giant and 
continuous accelerator (even if ultimately powered by super- 
nova explosions); the fact that Wommer et al. (2008) found 
they needed ^ 50 point, TeV 7-ray sources to reproduce the 
overall rather smooth distribution of 7-ray intensity over the 
GC region is consistent with such a picture. 

Over longer timescales, we have already listed above 
(§2.1.3) the evidence that the star-formation/supernova pro- 
cess itself is in steady state in the GC, certainly over the 
~30 Myr we estimate (following the work of Mo et al. 2010) 
necessary to ensure a steady supernova rate is reached for the 
lowest mass [Mzams = 8M0) SN progenitors and probably 
for much longer. 

3.2 Extensions to previous modelling 

We have significantly extended our treatment of the non- 
thermal phenomenology of the GC region to properly in- 
corporate other data covering the region and to introduce 
various physically-motivated constraints. Most importantly 
we explicitly require energy and mass conservation. Other 
constraints/data are explained below. One other extension to 
our previous modelling is that we add an extra term into our 
function related to the diffuse thermal X-ray luminosity 
of the region (see §3.8). 

3.3 Mass flows 

As described above, the SFR in the GC is constant when aver- 
aged over sufficiently long periods of time. Mass conservation 
then implies 
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MiN = MouT + SFR + McMz (1) 

where, given the relationship between SFR and the total 
surface mass density of gas (Kennicutt 1998), the total gas 
mass of the system is invariant on average, Mcmz = 0. M/jv 
represents all sources of mass falling on to the inner ~200 
pc (in diameter) of the Galaxy. This receives a substantial 
contribution from mass accreted through the plane but other 
sources of mass - including gas falling in (or back) from the 
halo above the GC and donated by the bulge's old stellar 
population - may also be important. In steady state, the 
difference between Mim and the SFR is the mass that must 
be expelled from the system on an outflow, Mqut- 

Note that it is the rate at which stars are formed - not 
the total number of stars - that is in steady-state. Never- 
theless, given that stars have finite lifetimes, the numbers of 
stars with Mzams sufficiently large that their lifetimes are 
less than the age of the system will be in steady state. Given 
the indications adduced above that the GC system has been 
forming stars at an approximately constant rate over a multi- 
Gyr timescale, the number of stars sufficiently massive to 
generate a supernova at death Mzams > SMq{t [8Mq] ~ 
30 Myr; e.g.. Mo et al. 2010) is in approximate steady state. 
In fact, the number of stars down to Mzams — 1 Mq is 
approximately constant (r [1 Mq] ~ 10 Gyr: Mo et al. 2010) 
given the system has been operating as long as we think it 
has. Of course, the number of lower mass stellar and sub- 
stellar objects (if they are formed) continues to grow with 
time as does the mass in compact stellar remnants. 

3.4 Power 

We assume that the mechanical power delivered by core 
collapse supernovae 

EouT = MsF X ( / M il}\M] dM J 

j-M^-p 

X / Esn[M] i/;[M] dM (2) 

(where tp = dN/dM is the IMF) drives the entire system with 
a (generally) subdominant contribution from stellar winds (of 
~ 10^^ erg/s; see Appendix B4 of Crocker et al. 2011b) at the 
10% level. Here the integrals are over, respectively the zero- 
age main sequence mass of progenitors, Mzams from the 
minimum mass object (quite possibly sub-stellar) into which 
cooling gas fragments, Mfrag, and from the minimum mass 
(at zero age) necessary for a star to explode as a core collapse 
supernova, Mdown, to the largest stellar mass arising from 
the star-formation process. Here we will assume that i) the 
fragmentation mass satisfies 0.07 ^ Mjrag/ Mq ^ 1.2 where 
the upper limit to M fmg is derived from the determination 
(Figer et al. 2004) that stars down to at least this Mzams 
exist in the GC; ii) the lower mass limit for a star to explode 
as a core collapse SN is ~ 8Mq (e.g., Smartt 2009); and iii) 
that the limiting upper stellar mass is 150 Mq as derived by 
Figer (2005). Note that our calculations are not very sensitive 
to the precise value of AIup- 

For simplicity in our modelling we ignore the sub- 
dominant power input from thermonuclear supernovae. There 
are other potential sources of power injection into the sys- 
tem (see appendix B6 of Crocker et al. 2011b, and references 



therein) which, however, are not particularly well constrained; 
amongst these we count processes associated with the super- 
massive black hole. As previously stated, one of the purposes 
of this paper is to demonstrate how the GC system can be 
kept ticking over without appealing to such processes. Note 
also that we expect for the GC environment that any radia- 
tion driving of the global outflow is negligible (cf. Thompson 
2009). 

3.5 Supernova energetics 

We need a prescription in our modelling for how the mechan- 
ical energy delivered by each supernova, Esn, evolves as a 
function of the zero-age mass of (Mzams) of the progenitor: 
here we will explore two limiting cases: 

(i) The 'standard' assumption that supernovae deliver a 
mechanical energy of 10^^ erg invariant with respect to the 
zero-age mass of (Mzams) of the progenitor. 

(ii) On the other hand, given there does seem to be some 
evidence for growth of Esn with Mzams (e.g., Nomoto et 
al. 2006, 2010; Utrobin & Chugai 2011) we wiU also explore 
a parameterization of this apparent growth (shown in fig. 2). 
We assume that the mechanical energy asymptotes to 6 x 
10^^ erg; this is simply set by the most energetic event (SN 
20031w) listed by Nomoto et al. (2010) rather than being, 
necessarily, some fundamental physical limit. Here (in the 
spirit of exploring the upper limit to the energy evolution) we 
ignore the possibility that the supernovae of some high-mass 
progenitors fizzle because of fall-back on to a newly-formed 
black hole (e.g., Nomoto et al. 2006) and we also ignore 
the fact that the mechanical energies arrived at by Nomoto 
et al. (2006) have been inferred under the assumption of 
spherical symmetry; the difference between the isotropic and 
real energies of any real SN is expected, however, to be less 
than a factor ~2 (see, e.g., supplementary material for Maeda 
et al. 2008). With the energy evolution parameterization we 
have assumed a hypernova (Esn 10^^ erg) requires a 
progenitor with Mzams > 26.5 M©. 



3.6 Power leaving system 

We require that the power being lost from the system, Equt, 
is matched by the power being fed into the system by core- 
collapse supernovae (and stellar winds), Ejn, 

EoUT = EjN 

EiN = Esn + Esw 

EoUT = Ethlrm + E^in + Eb + EcRp 

+ EcRe + E^f,f + Erad (3) 

Here Equt is composed of the following terms (cf. fig. 7): 
the i) thermal, Ethlrm — 4.1 fc_g T Mhot/rrip, and ii) kinetic 
power of the plasma in the outflow, E'^°^ = 0.5 Mhot Vhotj i") 
the magnetic energy advected by the outflow per second due 
to field lines frozen into the plasma, Eb — Eg^' Vhot/h 
where Eg^ is the total magnetic energy content of the 
GC region and /i ~ 42 pc its height; the power being 
lost into freshly-accelerated, non-thermal iv) protons and 
v) electrons; vi) the kinetic power of the cold, entrained 
gas in the outfiow, E^f^ = 0.5 Mcoid wLiti; s^nd vii) the 
power radiated by the plasma while it remains in the region. 
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Figure 2. A plot of the parameterization of the mechanical power 
released by a supernova as a function of the zero-age main sequence 
mass of the progenitor, Egpf[MzAMs]- The function is a rough 
parameterization of the data presented by Nomoto et al. (2010), 
disregarding the 'faint SN branch data points, and Utrobin & 
Chugai (2011) on SNIIP events. Note that parameterization of 
SN mechanical energy tops out at ~ 6 X 10^^ erg. Also notice 
that, under the influence of the Utrobin & Chugai (2011) data, 
at low MzAMS values the parameterized mechanical energy is 
considerably smaller than the canonical 10^^ erg. As explained in 
the text, the parameterization is effectively an upper limit to the 
real (population-averaged) dependence of Esn on Mzams- 



Erad = A[r, ne,n_r/] Vac where we adopt the parameteri- 
zation of Raymond et al. (1976) of the plasma volumetric 
cooling rate, A, as a function of temperature, electron and 
proton density. 



and mechanical power also being fed into non-thermal ISM 
phases). 



3.8 X-ray Luminosity of System 

We model energy and mass flows through the system and 
our modelling therefore gives us a handle on the steady-state 
plasma conditions in the GC which we can compare against 
data. We do not fix or constrain the temperature of the 
plasma beyond the empirically-motivated requirement that 
it be below 10* K, but we do calculate its 2-10 keV X-ray 
luminosity via both continuum free-free and line emission 
according to our parameterization of the results of Raymond 
et al. (1976). We add a term into our function accounting 
for the requirement that the predicted 2-10 keV luminosity 
reproduce the observed (Belmont et al. 2005) 4 x 10^^ erg/s 
(to within a factor 2 assumed to represent the la error in 
this luminosity). 



3.9 Number of Wolf-Rayet Stars in System 

The measured number of WR stars in the GC is 92 (Mauer- 
han et al. 2010; Hufimann et al. 2011). Allowing for the 
possibility that some such stars remain to be discovered, we 
conservatively assume that the steady-state number of WR 
stars predicted by our modelling should fall between this 
number and twice this number. The steady-state number of 
WR stars is given by 



Eoi 



MsF X 



M ililM] dM 



3.7 Velocity of Outflow 

Analytical solutions exist for determining the speed of a wind 
out of a starburst region as a function of distance given cer- 
tain simplifying assumptions (Chevalier 1992; Zirakashvili & 
Volk 2006; Strickland & Heckman 2009; Veilleux et al. 2005) 
like absence of halo drag, radiative losses, etc. We take a nu- 
merical approach as it best suits our purposes, in particular, 
allowing for the self-consistent treatment of non-thermal ISM 
phases (magnetic field, cosmic rays) which may themselves 
represent a significant sink of injected mechanical energy. 
Our approach allows the modelled system to distribute en- 
ergy between different phases without theoretical prejudice 
as to thermalization efficiency, cosmic ray acceleration effi- 
ciency, etc. Our modelling self-consistently determines both 
the mass loading and amount of cold gas entrainment that 
occurs in the outflow (these appear as free parameters in 
the analytical wind solutions; see §3.11 for the distinction 
between these). The outflow speed is given implicitly within 
our model by coupled equations, in particular, through the 
fact that the timescale over which the various ISM compo- 
nents are advected from the GC is one parameter controlling 
the power leaving the system (and the system is constrained 
to conserve energy globally) and also that this timescale is 
also a controlling parameter for the steady-state non-thermal 
particle spectra (and, therefore, their consequent radiation). 
Note that, in general, we flnd below that - as it leaves the 
boundary of the GC system - the outflow is significantly 
sub-sonic (as a result of both mass loading and entrainment 



i^[M] twr[M] dM 



(4) 



where twr[Mzaais] is the duration of the WR phase experi- 
enced by a star of given Mzams and M^^^ is the minimum 
mass required in order that a star experience such a phase. 
Here there are considerable theoretical uncertainties influ- 
enced by whether, at the population level, the most important 
channel for lifetime-integrated massive stellar mass loss is via 
(single star) winds, eruptions, or binary mass transfer. This 
is very much a held of active research and debate (Dessart 
et al. 2011; Modjaz et al. 2011; Meynet et al. 2011; Smith et 
al. 2011; Boissier & Prantzos 2009; Smartt 2009; Eldridge et 
al. 2008; Meynet & Maeder 2003); for deflniteness we shall 
adopt a parameterization of the results of Meynet & Maeder 
(2003) for the massive, single-star evolution of non-rotating 
and stars rotating initially at 300 km/s. Rotating stars lose 
mass more quickly and get to spend longer in the WR phase 
before core collapse, so evolution with rotation predicts a 
larger steady-state number of WR stars for a given SFR 
and IMF than the no rotation case are require that our 
modelling predict. Within the single star evolution paradigm 
M^'^^ ~ 25 Mq . To bracket the theoretical uncertainties we 
require that NwR{no rot) < 2 x 92 and NwR{rot) > 92. 

3.10 Other constraints 

Other constraints we enforce are as follows: 

Pree-free emission from system: Our detailed modelling 
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accounts for thermal bremsstrahlung (free-free) emission 
from individual Hll regions and/or dispersed plasma within 
the GC in fitting to the higher-frequency part of the radio 
continuum spectrum (and self-consistently accounts for free- 
free absorption of radiation at lower frequencies). Free- free 
emission can be used to derive an empirically-calibrated lower 
limit (because of dust-absorption) to the SFR (see Eq. 11 of 
Murphy et al. 2011) in the system and we demand that the 
SFR and free-free satisfy this implicit constraint. 
Total infrared emission from system: Ltir offers a 
rather robust, empirically-calibrated (but IMF-dependent) 
handle on the SFR in the system. We do not directly model 
the system's total infrared output but we do conservatively 
require that the modelled SFR (of > 5 Mq stars) predict 
(Condon 1992) a Ltir no more than a factor of two larger 
than the output observed (1.6 x 10*^ erg/s Launhardt et al. 
2002) from the system. 

Ejecta masses: We require that the masses ejected by the 
SN explosions ending stars' lives (accounting for modelled 
mass loss over their lifetimes and conservatively assuming 
a 1.4 Mq compact remnant irrespective of Mzams) be 
positive. 

Explaining radio continuum from GCL: As introduced 
above, there is good evidence that non-thermal electrons 
injected at the plane are carried out of the immediate GC re- 
gion to synchrotron-illuminate the GCL. In fact, as discussed, 
radio observations (LaRosa et al. 2005; Crocker et al. 2010a) 
reveal a distinct, extended, non-thermal radio continuum 
source on even larger scales around the GC (6° in longitudi- 
nal extent) of which the GCL non-thermal emission forms 
only a part. Our previous modelling has shown that the inner 
200 pc region can supply enough power (~ 10"^* erg/s) in 
hard-spectrum cosmic ray electrons to explain these obser- 
vations (and this is borne-out by the modelling presented 
below); the only question is whether the outflow transports 
this population on timescales short enough with respect to 
the electrons' loss time in the dense and highly-magnetized 
GC ISM. Given the physical plausibility of the requirement 
that the GC supply the GCL electrons, we therefore demand 
that the modelled wind be able to transport the advected elec- 
tron population to 100 pc over the electrons' loss time (at 
the energy corresponding to synchrotron-emission - for the 
modelled magnetic field - into the highest frequency range 
where emission is still manifestly non-thermal in character, 
viz. 10 GHz). 

3.11 Collimation, mass loading, and entrainment 

Following Strickland & Heckman (2009), we distinguish here 
and in our modelling between (centralised) mass loading and 
mass entrainment. In the former, additional ISM material 
is heated to plasma temperatures and co-mixed with the 
supernova and stellar wind material in the energetically 
dominant 'wind fluid' (Strickland & Stevens 2000) thereby 
affecting the entire flow in a global, distributed process. In the 
latter, local process, a cooler and denser gas phase, stripped 
from ambient cold gas by the ram pressure of the wind fluid, 
is carried along with the flow but remains a distinct phase. 
Note the entrained gas may actually come to dominate the 
mass efflux. Physically, of course, a real wind or outflow is 
characterised by multiple co-evolving phases and entrainment, 
e.g., may lead to loading but we ignore these subtleties in 



our modelling. The physical picture suggested to us by the 
data and which we advocate is that the dense and massive, 
star-forming gas ring (Molinari et al. 2011) coUimates the 
wind outflow (cf. Zubovas et al. 2011; Strickland & Stevens 
2000; Westmoquette et al. 2011; the inferred large-scale, 
poloidal field structure may also have a role in collimating 
the outflow; Morris 2007). Indeed, the radio continuum data 
suggest that the outflow has a projected width matching 
the projected radius of this ring (see flg.l). Conversely, ram 
pressure stripping of H2 from the inner edge of the ring by 
the escaping plasma outflow supplies gas to be loaded on to 
the outflow. 

We quantify gas entrainment following the prescription 
of Martin (2005) who shows that the ram pressure of the 
hot outflow, of density pwind and speed Vwind, accelerates 
entrained, cold gas clouds, of density pc, to a terminal speed 
given by 



We assume in our modelling that the density of the cold 
material is that corresponding to the number density which 
the non-thermal particle populations are primarily sampling 
and within which they generate the non-thermal radiation 
we detect. The entrained gas can be expected to be shock 
heated by the plasma wind fluid, some fraction potentially 
to X-ray emitting temperatures (e.g., Strickland et al. 1997) 
and potentially converting from 'entrained' to 'loaded'; we 
leave a detailed treatment of this for future work.. 



4 RESULTS 

Given the above model and constraints we delimit the param- 
eter space providing a good description of the GC environ- 
ment using minimisation to fit to the diffuse, broadband 
emission detected from the GC. We are fit to 6 radio data 
points (at 74 MHz, 1.4 GHz, 2.4 GHz, 2.7 GHz, 8.4 GHz, 
and 10.3 GHz), 9 HESS 7-ray points (from 2.7 x lO" eV 
to 1.3 X 10" eV), and one diffuse X-ray flux datum (2 to 
10 keV); there are 9 fltting parameters and one constraint 
(energy conservation) so 8 degrees of freedom overall. We 
show an example fitted broadband spectrum in Appendix B. 

As discussed, we employ two prescriptions for the me- 
chanical power delivered by core-collapse supernovae into the 
ISM to bracket the reasonable possibilities, viz.: that this 
be invariant at 10^^ erg per SN and that this be growing 
function of Mzams as described in §3.5. 

In general, we find fits acceptable at 2a confidence for 
the control parameter Mjn in the ranges 0.1 — 2MQ/year 
and 0.2 — 5M0/year for these two cases (fig. 3). 



4.1 Star formation rate 

On the basis of our constrained modelling, we determine a 
SFR in the system that, at 2cr-level, lies in the range 0.04- 
0.12 Mq /year (fig. 4) with best-fit values of 0.08 and 0.12 
Mo/year for the cases, respectively, of Esn = 10^^ and 
Esn[Mzams]- These values agree well with the previous, 
independent determinations set out in §2.1.1. 
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Figure 3. ^ ^ function of the mass flux into the GC region, 
Mjff, for the case of (sohd,blue) Eg^ = 10^^ erg (constant with 
respect to Mzams) a-nd (dashed, purple) Esn an increasing 
function of Mzams as described in §3.5. 
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Figure 4. Star formation rate as a function of the mass flux into 
the GC region, M/jv, for the case of (solid, blue) Esn = 10^^ erg 
(constant with respect to Mzams) and (dashed, purple) Eg^ 
an increasing function of Mzams as described in §3.5. In this 
figure and the following we show modelled parameters only over 
the range of M/jv where they are acceptable at 2cr-level. 

4-1.1 Inferred age of system 

Fig. 5 displays the time required by the system to form the 
~ 8 X 10® Mq stellar population (Launhardt et al. 2002) of 
the R < 120 pc region of the Nuclear Bulge being modelled. 
Given an age to the system of ^ 10 Gyr, this indicates 
minimum values for Min of ~ 0.4 and 2MQ/year for the 
cases, respectively, of Esn = 10®^ and Esn[Mzams]- Given 
we also find upper limits (at the 2a level) for the SFR 
(corresponding to the largest plotted values of Mjn for both 
Msf curves) corresponding to a formation time of ~6 Gyr, 
this latter constitutes the minimum age of the system within 
our scenario. 

4.2 Modelled Powrer 

Our modelling (fig. 6) reveals that the inner 200 pc region 
requires an input power of ^ 2 x lO*"^ erg/s to sustain its non- 
thermal radiation and, more importantly energetically, the 
inferred outfiows of non-thermal and thermal ISM phases. 
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Figure 5. The total time required for the GC to supply the 
~ 8 X 10* Mq stellar population (Launhardt ct al. 2002) of the 
R < 120 pc region of the Nuclear Bulge for the case of (solid, blue) 
Esn = 10^"^ erg (constant with respect to Mzams) and (dashed, 
purple) Esn an increasing function of Mzams as described in 
§3.5. The horizontal, dotted, yellow line shows the nominal, 10 
Gyr age of the system. 
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Figure 6. Ejn- Total power being delivered into the GC by 
supernovae and, sub-dominantly, stellar winds for the case of 
(solid, blue) Esn = 10^^ erg (constant with respect to Mzams) 
and (dashed, purple) Esn an increasing function of Mzams as 
described in §3.5. 

One important implication of this is that the finding by 
Yasuda ct al. (2008) that the current level of star formation 
in the GC is low (with the total FIR being dominated by old 
K and M giants) seems difficult to sustain (i.e., our results 
are consistent with the SFR inferred were young, massive 
stars to dominate the system's radiative output as is the 
case for star-burst-like environments and more generally; e.g. 
Thompson ct al. 2006). 

Our modelling also allows us to determine the efficiency 
with which the mechanical energy injected by supernovae 
is, on the one hand, converted into heating or moving the 
ISM (their 'thermalization efficiency') and, on the other, 
converted into freshly-accelerated, non-thermal particles. In 
common with some analyses of the dynamics of winds driven 
out of star-forming nuclear regions of external galaxies (e.g., 
Strickland & Heckman 2009) we find a high thermalization 
efficiency for the GC system for both the Esn = 10^^ erg 
and Esn[Mzams] cases. 



14 Roland M. Crocker et al. 





40 




39 


erg/s 


38 






t 


37 


-1 


36 




35 




34 



j^*^.'-T *** : 



-kinetic hot 
kinetic cool 
- thermal 

■ magnetic 
■radiative 

■ P 



-1.0 -0.8 -0.6 -04 -0.2 0.0 0.2 
Log[M,N/M^/yr] 

Figure 7. Equt- Power being delivered into various phases for 
the case of Egjv = 10^^ erg (constant with respect to Mzams) 
as a function of the mass flux into the GC region, Min- The 
curves, as denoted in the legend, are for i) the kinetic power of 
the hot plasma in the outflow, ii) the kinetic power of the cool gas 
entrained by the outflow, iii) the thermal power advected by the 
hot plasma in the outflow, iv) the magnetic energy advected by 
the outflow per second due to field lines frozen into the plasma, 

v) plasma radiative losses in the GC region, and power being 
delivered into the freshly-accelerated populations of non-thermal 

vi) protons and vii) electrons Note that the power going into 
the freshly-accelerated proton population is almost invariant at 
~ lO^S erg/s over the well-fitting region of M/]v; this is precisely 
the proton power required, in steady state, to sustain the 7-ray 
luminosity of the Fermi Bubbles. Also note that the power going 
into the freshly-accelerated electron population is almost invariant 
at ~ 10'^* erg/s over the well-fitting region of Min\ this is the 
electron power required (Crocker et al. 2010a), in steady state, to 
sustain the synchrotron luminosity of the extended, non-thermal 
emission detected (LaRosa et al. 2005) around the GC including 
the GCL structure. 
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Figure 8. Equt- Power being delivered into various phases for 
the case of Esff[MzAMs]- The curves are as explained in the 
caption to fig. 7 



-0.2 

>; -0.4 

I ^ 

S -0.8 
o 

-1.0 
-1.2 




■ thermalization eff . 

■ CR eff. 



-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 
Log[M,N/M /yr] 



Figure 9. Modelled thermalization efficiency and cosmic ray 
acceleration for the case Esn = lO^"*^ erg. 
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Figure 10. Modelled thermalization efficiency and cosmic ray 
acceleration efficiency for the case Egif[MzAMs]- 



With respect to the efficiency of the system as a cosmic 
ray accelerator, we can say that, were the Esn = 10^^ erg 
prescription correct, the system would be typically or close 
to typically (e.g., Hillas 2005) efficient, losing 5-10% of power 
into these non-thermal particles. Were the Esn[Mzams] 
prescription correct, a rather low cosmic ray acceleration 
efficiency of ^5: 1% would be implied. 



4-2.1 Advected magnetic field 

Figs. 7 and 8 make clear that magnetic energy losses in the 
form of field-lines frozen into the advected plasma contribute 
substantially to the total energy budget of the system^. The 
magnetic energy losses are, in fact, larger than or comparable 
to the kinetic power of the outflow indicating that the mag- 
netic field is an important determinant in the dynamics of the 
outflow (e.g. Beck et al. 1996). This is consistent with, on ~ 
degree scales, the complex but ordered phenomenology of the 
outflowing dust filaments revealed by MIR observations (as 
briefly reviewed above). On tens of degree scales, magnetic 
field must equally play an important part in governing the 
evolution of the Fermi Bubbles as revealed by the detection 
of polarised microwave radiation from the GCS and part of 
the edge of the Northern Bubble (Jones et al. 2012). We will 
consider the evolution of the magnetic held - governed by 
reconnection/relaxation and adiabatic losses - between the 
scales of the GC region and the full size of the Bubbles more 
fully below. 



4.3 Mass Outflows 

A first important point here is that, although the entrained, 
cooler gas is moving more slowly than the hot outflow, it 
represents a larger mass flux (fig. 16). We deal with these 
different phases separately below. 



^ In principle, the calculated energy loss into the magnetic field 
disregards the fact that the vertical outfiow will not advect the 
poloidal magnetic field component and it therefore represents a 
de facto upper limit on the true losses. Giiven, however, that the 
outfiow is entraining cooler gas and frozen-in field lines and - for 
the dense gas — the field lines have been sheared into a toroidal 
configuration consistent with polarimetry measurements (Chuss 
et al. 2003) we expect the real magnetic energy loss rate to be 
close to that displayed. We thank Mark Morris for bringing up 
this point. 
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Figure 11. Speed of hot (solid) and cold (dashed) outflows for 
(blue) the case of Esn = 10^^ erg; and (purple) Esn an increasing 
function of Mzams described in §3.5 
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Figure 12. Mass flux in hot (solid) and cold (dashed) outflows 
for (blue) the case of Eg^f = 10^^ erg; and (purple) i?sjv an 
increasing function of Mzams as described in §3.5 

4.3.1 Hot Outflow 

We calculate the height that the escaping hot material might 
reach under the assumption that it is moving ballistically 
after crossing the boundary of the region. We assume that 
the strong gravitational potential of the region is as given by 
the parameterization of Breitschwerdt et al. (1991) for the 
case of vanishing Galactocentric radius. This gives an escape 
speed of ~ 810 km/s for material launched from z; = 42 pc 
(the half-height of our region) . 

Given the strong gravitational deceleration, ballistically- 
moving hot material only reaches heights of 300 - 1000 
pc for the case of Esn = const, with considerably higher 
ranges predicted for the Esn[Mzams\ case. This, however, is 
likely to be an underestimate of the true height the material 
reaches: the pressure gradient represented by thermal and 
non-thermal components will likely further accelerate the 
outflow after it passes the region's boundary. Here a rough 
upper limit - accounting for this effect - is given by the 
assumption that the material moves at the sound speed: 
this would imply that the outflow is not gravitationally 
bound for Min < 0.4 /year and 3 /year for the cases, 
respectively, of Esn = const and Esn[Mzams]', however, 
these low ranges of Mjn are disfavoured by the consideration 
that they result in SFRs too small to fill-out the Nuclear 
Bulge stellar population over the age of the Galaxy as already 
discussed (cf. fig. 5). Note that for the favoured Mjn range 
and the Esn = 10^^ erg case the outflow velocity at the 
region's boundary and the sound speed differ by less than a 
factor of 3. 

A rather firm conclusion of our modelling, then, is that 
the material leaving the GC does not escape to infinity, i.e., 
the outflow is not a true wind but rather a fountain (Breg- 
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Figure 13. The speed of the hot GC outflow for the case of 
(solid, blue) i?sjv = 10^^ erg (constant with respect to Mzams) 
and (dashed, purple) Esn an increasing function of Mzams 
as described in §3.5. The dotted, yellow horizontal line shows 
the escape velocity for material launched from a height of 42 pc 
given our assumed gravitational potential. The blue and purple 
dot-dashed lines show the sound speeds in the plasma given its 
modelled temperature for, respectively, the Esn = 10^^ erg and 
EsnIMzams] cases. The long dashed curves likewise show the 
Aflven speed for the same two cases. Note that the outflow has only 
just reached its Alvenic point on the boundary of the modelled 
region for favoured values of M/jv; inside this region, the magnetic 
field will rotate rigidly, consistent with the phenomenology of the 
region's non-thermal filaments (Yusef-Zadeh et al. 1987; Morris & 
Yuscf-Zadch 1989). 
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Figure 14. The height attained by the hot material in the outflow 
assuming it is moving ballistically after crossing the boundary of 
the region. This is unlikely to be true: the pressure represented by 
thermal and non-thermal components will likely further accelerate 
the outflow and it will therefore reach greater heights than plotted, 
which therefore represent lower limits. A rough upper limit to the 
height of the outflow is presented for the case that the outflowing 
material is assumed to cross the boundary at the sound speed 
(in our assumed gravitational potential the escape speed is ~ 810 
km/s at 42 pc height; if the upper limit on the outflow height is not 
plotted the sound speed is above 810 km/s). The case of Esn = 
10^^ erg (constant with respect to Mzams) is shown in blue for 
(solid) the minimum height and (dot-dashed) the maximum height; 
the case of Esn an increasing function of Mzams as described 
in §3.5 is shown in purple for (dashed) the minimum height and 
(dot-dashed) the maximum height. 
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Figure 15. Inferred amount of centralised mass loading for the 
case of (solid, blue) i?siv = 10^^ erg (constant with respect to 
MzAMs) s-iid (dashed, purple) iJgjv a-n increasing function of 
l^ZAMS a-s described in §3.5. 



4-3.2 Centralised mass loading 

We can calculate the amount of centralised mass-loading, 
(3, within our model (fig. 15); this is the ratio between the 
total hot mass efflux (which includes swept-up, heated ISM 
gas) and the directly injected hot gas originating as super- 
nova ejecta and stellar winds. We find /3 ~ 3-10 for the 
EsN = W^^ erg case over the favoured region of parame- 
ter space (log [M/iv/M© /year] > - 0.4) and /3 ^ 10 for the 
case of Esn[Mzams]- In comparison, the recent study of 
Strickland & Heckman (2009) finds that /3 lies in the range 
1-3 for the star-burst conditions in M82, with a practical 
upper limit, for this system at least, at ~ 10. Other studies 
have countenanced or suggested somewhat higher values for 
star-burst environments (e.g., ~ 5 according to Suchkov et 
al. 1996, ~ 10 according to Martin 2005). Taking these inde- 
pendent estimates at face value, the Esn[Mzams] scenario 
does appear to be disfavoured by the large amount of mass 
loading it requires. 



man 1980). Qualitatively, this matches the UV-absorption 
phenomenology introduced above (cf. §2.3.6) which demon- 
strates the existence of highly-ionized material fountaining 
up to heights of 12 ± 1 kpc (Keeney et al. 2006), a scale 
interesting close to the ~ 10 kpc heights of the Fermi Bub- 
bles. The launching of material to such heights requires a 
mechanism distinct from - and significantly more powerful 
than - the 'standard' disk-halo gas connection mediated by 
blow-out of super bubbles around disk supernova associations 
(de Gouveia Dal Pino et al. 2010) and even seems difficult 
to explain within our model for the Esn = 10^^ erg case 
(as fig. 14 demonstrates). If we demand that the global out- 
fiow shoot material directly to ~ 10 kpc then it would seem 
Esn [Mzams] is preferred. 

However, even if the global ' super '-outfiow we model is 
incapable of reaching these sort of distances directly, there 
are two effects that may mean GC material still reaches these 
heights: 

(i) The activity of individual mini star-bursting events - 
leading to the creation of the GC's super stellar clusters (of 
which the GC, Arches, and Quintuplet are merely the most 
recent examples) - might achieve this (Rodriguez-Gonzalez 
et al. 2009). The GC Spur may represent evidence for just 
such an outfiow-within-an-outflow as previously discussed 
(§2.3.7). 

(ii) The ejected low-density thermal and CR plasma may 
rise buoyantly once ejected into the bulge (Su et al. 2010, 
also cf. Rodriguez-Gonzalez et al. 2009). X-ray and radio 
observations (e.g., McNamara et al. 2005) reveal that such 
a mechanism certainly operates on the scale of galaxy clus- 
ters; the rather slow velocities involved would, moreover, be 
entirely consistent with the general expectation afforded by 
our previous work (Crocker & Aharonian 2011) that the 
formation timescale for the Fermi Bubbles is rather long. In 
this regard, we note the work of Raley et al. (2007) who 
modelled the buoyant ascent of bubbles formed by individual 
supernovae in the Galactic halo and found surprising slow 
rise speeds of only 5-15 km/s (not too much larger than 1-2 
km/s Bubble growth speed implicit in our scenario) taking 
into account a very high, effective drag coefficient for the 
bubbles in the halo plasma. 



4.3.3 Cold Outflow 

We find a mass flux in the cold outflow 2-6 times larger 
than in the hot outflow over the favoured Min parameter 
space (flg 16). The entrained cold gas is (relatively) slow 
moving and only reaches heights of a few x 100 pc from the 
plane (at most). It will, therefore, fall back on the GC. Thus, 
while the fountaining cold gas seems to dominate the mass 
flux out of the system, it also makes a large contribution 
to the total mass flux into the system, (Min) and the net 
mass fiux represented by the cold gas is small. Still, this 
cold gas circulation is likely to have important dynamical 
effects. Firstly, the fall of this material back on to the plane 
likely represent a significant source of turbulent stirring of 
the gas there (e.g. Klessen & Hennebelle 2010), additional 
to the in-situ effect of expanding Hll bubbles, stellar winds, 
and supernovae consequent to the region's star formation^". 
Secondly - as we describe below - the ejection of cold, high- 
metallicity gas into the halo and its subsequent mixing with 
the in-situ plasma should 'catalyse' the further condensation 
and accretion of plasma out of the halo (Marinacci et al. 
2010, 2011; Binney & Fraternah 2011). 

4.3.4 Cold gas filling factor 

The modelled filling factor of the entrained cold gas is 
shown in fig. 20. For the case of both Esn = const and 
Esn[Mzams] the filling factor is rather large over some of 
the favoured range of the control parameter Min with re- 
spect to the expectation from star-burst winds (which would 
suggest a range 0.1 - 0.01; e.g., Strickland et al. 1997). This 
may indicate a break-down of the model assumption that 
the non-thermal emission dominantly arise in the (relatively) 
cold, entrained material. Alternatively, the filling factor re- 
sult may be correct and indicate a point of difference between 
the GC situation and true star-burst environments. 

The kinetic power in the entrained gas approaches 10"^^ erg/s 
for the favoured Mjn range — cf. fig 7. We previously estimated 
the power lost into turbulence dissipation in the region at 4 X 10^® 
erg/s: see Appendix B 2.5 of Crocker et al. (2011b). 
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Figure 16. The mass flux in the cold (entrained) outflow relative 
to the hot outflow. 
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Figure 17. The height attained by entrained, cold gas after being 
launched from 42 pc height were it to move a ballistic trajectory 
for the case of (solid, blue) Egfj = 10^^ erg (constant with respect 
to MzAMs) ^iid (dashed, purple) Egj^ an increasing function of 
f^ZAMS described in §3.5.. Recent Nanten-II measurements 
(see §2.3.3) indicate a total mass of few X 10® Mq of molecular gas 
in an extended (100-200 pc) halo around the GC region. 

One piece of evidence tending one to the latter view is 
that analysis of the large columns observed in absorption 
towards the region (Goto et al. 2008; Yusef-Zadeh et al. 2007; 
Oka et al. 2005) and different emission lines of i^CO (Oka 
et al. 1998) and other molecules (Rodn'gucz-Fcrnandcz et al. 
2001) indicates the presence of a highly-ionized ( ^ 10~^), 
comparatively diffuse (~ 100 cm~^) and hot (~ 250 K) 
molecular phase. This 'envelope' H2, which represents ~30% 
of the total molecular gas by mass (Ferriere et al. 2007), 
appears to be unique within the Galaxy and has been claimed 
to have a high filling factor, perhaps approaching 100 % (Goto 
et al. 2008). 

4.4 Inferred Mass Accretion Rate 

As revealed above, a number of indicators come together to 
suggest that our control parameter, M/jv - the total mass 
being fed into the system - has a lower limit at around 0.4 
M0/year: i) given the star-formation rate determined by our 
modelling as function of Mjjv, we find that to fill out the 
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Figure 18. The steady state mass of cold gas which is either 
entrained by the outflow or falling back to the plane after reaching 
its maximum height for the case of (solid, blue) Esn = 10^^ erg 
(constant with respect to Mzams) ^iid (dashed, purple) i?sjv 
an increasing function of Mzams a-s described in §3.5. Recent 
Nanten-II measurements indicate a total mass of few X 10® Mq of 
molecular gas in an extended (100-200 pc) halo around the GC 
region (see §2.3.3); this represents independent support to the 
notion that Af/jv ~ 1 Mq /yeai. 
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Figure 19. Modelled density of the phase wherein the observed 
non-thermal radiation is generated, assumed to be identical to the 
cold gas entrained by the outflow. 



stellar population of the inner regions of the nuclear bulge 
in ~ 10 Gyr at least this Mjn is required; ii) to explain the 
amount of mass in the GC lobe at least this Mjn is required, 
more generally, the recently-identified ^ Mq molecular 
halo around the CMZ, seems to require Mjn ~ 0.4 /year; 
iii) if we require that the GC sit on the Kennicutt (1998)- 
Schmidt relation we seem to require Min ~ OA Mq /year; 
and iv) the requirement that the Bubbles be slightly pos- 
itively (or even neutrally) buoyant in our scenario implies 
r] ^ 1 which is again satisfied for Min ~ (0.6 — 1) Mo/year 
for the Esn = 10^^ erg case. 

We find a 2a upper limit on Min (for the Esn = 10^^ 
case) at 1.8 A/Q/year. It is intriguing that over the favoured 
Min range for the Esn = 10^^ case there are very similar 
energy loss rates into the kinetic power of the hot outflow, 
cosmic ray acceleration, and advected magnetic field (cf. fig. 7 
and 8). 
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Figure 20. Inferred filling factor of the cold gas entrained by 
the outflow (also assumed to be the phase wherein the observed 
non-thermal radiation is generated). 
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Figure 21. Mass flows for the case of Esjv = 10^^ erg (constant 
with respect to Mzams)- Given that the cold gas is likely to 
fountain back on to the system, we can determine an upper limit 
on the true, net inflow into the system M^^=*" < SFR + M^ot- 
Following the work of Marinacci et al. (2011) we also display 20% 
of the cold mass flow as a rough estimator of the scale of the 
net mass growth (with every cycling of the cold gas) due to the 
phenomenon of high-metalicity, outflowing Hi mixing with and 
subsequently 'catalysing' the condensation and accretion of halo 
plasma. 

Our modelling plus the requirement that the GC star- 
formation fill out the stellar population of the Nuclear Bulge 
over 10 Gyr or less point to a total gas mass that has been 
processed through the GC over this timescale of ^ (3 — 
10) X lO^'M©. The mass flux out of the system is dominated 
by the cold, entrained material but we have already seen 
that this material fountains back on to the GC. As this 
material is cycling - rather than being truly ejected - Min = 
SFR+MouT = SFR+Mhot+Mcoid would represent an over- 
estimate of the mass accretion rate required to sustain the 
system in steady state. This is, instead, given approximately 
by Mi^'^" SFR + Mhot- Note this is an still upper limit 
because - given that even the hot outflow does not reach 
the escape speed - some proportion of this plasma (given by 
the details of the interaction between the hot outflow and 
halo) will also fountain back on to the region. We plot the 
upper limit on Alf^*" in figs. 21 and 22 for the two cases 
of EsN investigated: it is ~ 0.5 /year over the favoured 
parameter space for both these cases. From fig. 23 we can 
infer that 3 x 10^ Mq of gas must be accreted on to the 
GC over 10 Gyr. We discuss below the implications of this. 
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Figure 22. Mass flows for the case of Egpf an increasing function 
of Mzams described in §3.5. 
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Figure 23. Total masses assembled/processed over an assumed 
10^*^ years and for the i?sjv = 10^"*^ erg case. 

4.5 GC ISM conditions 

4-5.1 Magnetic field 

Consistent with previous work (Crocker et al. 2010a, 2011a; 
Crocker & Aharonian 2011; also see Contini 2011; Spergel 
& Blitz 1992) we find our modelling favours a very high 
magnetic field amplitude in the GC, of order 100 /xG (i.e., 
~400 times the energy density of the local ISM field). As 
we discuss above and below, such a strong field can have 
important gas dynamical effects. Another very important 
aspect of GC magnetic phenomenology is the existence of the 
non-thermal filaments (NTFs). These are thin, synchrotron- 
illuminated structures running mostly perpendicular to the 
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Figure 24. Inferred magnetic field amplitude in the GC ISM for 

the cases of (solid, blue) Egpf = 10^^ erg (constant with respect 
to Mzams) ^-i^d (dashed, purple) Egpf an increasing function of 
Mzams described in §3.5. 
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Figure 25. Temperature of the diffuse plasma in the GC for the 
case of (sohdjblue) -Egjv = 10^^ erg (constant with respect to 
!^ZAMs) (dashed, purple) -Esjv an increasing function of 
^ZAMS a-s described in §3.5. 

Galactic disk and characterised by very high magnetic field 
amplitudes, ~ mG (Yusef-Zadeh et al. 1987; Morris & Yusef- 
Zadeh 1989). Given that we have claimed that it is enhanced 
star formation in the GC region that drives all the non- 
thermal phenomenology dealt with here and a large-scale 
outflow, it is intriguing that the NTFs are unique within the 
Galaxy to the longitude range of this region. In fact, magnetic 
field advection by the wind may well have a crucial role in 
forming the NTFs as analogues to cometary plasma tails, 
formed from the interaction of the large-scale, magnetized 
plasma outflow draping the region's dense molecular clouds. 
(Shore & LaRosa 1999; Boldyrev & Yusef-Zadeh 2006). 

4-5.2 Plasma conditions 

Our modelling indicates that it is possible for the system to 
heat the outflowing plasma to the 'very-hot' temperatures 
~ 7 X 10^ K pointed to by X-ray observations. However, 
that the plasma be this hot is clearly disfavoured for the 
EsN = 10^^ erg case at the sort of Min values suggested 
by the considerations given above; even the Esn[Mzams] 
case favours somewhat cooler temperatures over some of its 
favoured Min range. 



5 IMPLICATIONS 

5.1 Modes of Gas Accretion on to Inner Galaxy 

Chandran et al. (2000) estimate - to order of magnitude 
precision - that 3 x 10^ Mq of gas has fallen on the central ~ 
200 pc over the Galaxy's lifetime; this would imply an average 
accretion rate of ~ 0.3 Mq /year. Figer et al. (2004) estimate 
a mass accretion rate of ~ 0.4 /year from dividing the 
mass in molecular ring circumscribing the CMZ, 8 x 10® Mq, 
by its orbital period, 2 x 10^ year. A conservative upper limit 
on the accretion rate can be derived from the assumption that 
all the material falling on the Galactic disk, < 10 Mq /year 
(Combes 2004) - and more probably 0.5 — 5 M© /year (e.g. 
Wakker et al. 1999; Klessen & Hennebelle 2010, and references 
therein)- finds its way to the GC. 

Accretion through the disk: The overall position and 



size of the CMZ is presumably controlled by the gravita- 
tional dynamics dictated by the Galactic bar (Binney et 
al. 1991; Morris & Serabyn 1996; Stark et al. 2004). The 
non-axisymmetric gravitational potential of the bar induces 
torques on disk gas that cause it to fall inwards at a rate 
estimated (Morris & Serabyn 1996) at 0.1-1 MQ/year. 

Other mechanisms could, however, be acting to supple- 
ment accretion of gas on to the GC and may be particularly 
important in achieving further transfer of gas inward of the 
characteristic ~100 pc radius of the X2 orbit family. Such 
additional mechanisms include the action of a secondary, 
nested bar (Namckata et al. 2009), dynamical friction, and 
shocks associated with the XI — >■ X2 orbit transition (see 
Morris & Serabyn 1996, §3.2 for a more inclusive list). Also 
of interest are: 

Magnetic torquing: If sufficiently strong, the GC magnetic 
field might have appreciable effects on gas dynamics in the 
region (Morris & Serabyn 1996). In particular, magnetic 
viscosity may provide a channel for accretion of gas to radii 
smaller than that of the molecular ring established by the 
larger-scale gravitational dynamics. Scaling the results of 
Balbus & Hawley (1998) and Beck et al. (1999) (also see 
Beck et al. 2005), we find that magnetic stress should provide 
a mass infiow rate inside 100 pc of Mim™^ — h Q,~^ ~ 
0.2 Mq /year, where the total magnetic field amplitude is 
assumed to be _B = 100 /iG, the gas scale height ~13 pc 
(Fcrriere et al. 2007) and Q, ~ 100 km/s/100 pc is the angular 
velocity. This is interestingly comparable to other accretion 
channels. 

Bulge stellar mass loss: A minimum level of ~ 
O.2M0/year (Jungwiert et al. 2001) of gas is supplied to the 
outer bar (Stark et al. 2004) by mass loss from the evolved 
bulge stars. Accretion of material out of the slowly-rotating 
Galactic bulge represents a dilution of disk specific angular 
momentum (Morris & Serabyn 1996, and references therein). 

We have already seen (§2.3.4) that there is interesting 
evidence for a small amount of relatively pristine gas in the 
GC region (Lubowich et al. 2000), especially in the outer 
XI orbit family or in the process of being transferred from 
these outer orbits to the inner X2 orbits (Riquelme et al. 
2010; Jones et al. 2011). This represents rather compelling 
evidence for the accretion of material out of the halo on to 
the GC which would require a mechanism or mechanisms 
distinct from those listed above. 

Three, non-mutually-exclusive mechanisms might oper- 
ate to collect this gas, all potentially related to the outfiow we 
have identified here and, therefore, potentially self-catalysing. 

Dust cooling: Ambient dust grains provide for the colli- 
sional cooling of plasmas (e.g., Montier & Giard 2004; Natale 
et al. 2010). The operation of this mechanism in the halo 
above the GC requires the injection of dust into the region. 
Mid-infrared maps trace outfiows of dust coincident with 
radio continuum spurs running north and south from the GC 
region (Bland-Hawthorn & Cohen 2003; Stolovy et al. 2006; 
Morris et al. 2006) as already noted (§2.3.3), thus there is 
good evidence for entrainment of dust in GC outfiows and, 
presumably, its injection into the halo above the GC where 
it may serve to catalyse the cooling of the halo plasma; we 
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leave a quantitive treatment of this potential mechanism to 
further work^^. 

Shocks ahead of Giant Molecular Loops: Also as al- 
ready noted, observations with the Nanten-II instrument 
(Fukui et al. 2006) suggest the presence of giant, rising loops 
of molecular material in the GC region. The somewhat sur- 
prising presence of large amounts of molecular material arch- 
ing over the entire length of each of these loops (rather than 
being strongly concentrated at the 'foot-points' where the 
gas might be expected to fall subsequent to each loop's rise 
under magnetic flotation) may simply be a question of timing 
but it does seem, in general, to place the magnetic floatation 
via the Parker instability interpretation under some strain 
(Morris et al. 2006). Riquclme et al. (2010), moreover, have 
also found relatively nuclear-unprocessed molecular gas coin- 
cident with the GML footpoints, which again is somewhat 
surprising if the loops represent molecular gas magnetically- 
levitated out of the disk close to the GC (for, were this true, 
the footprint isotopologue ratios would be consistent with 
the highly-processed gas found in the rest of the inner CMZ). 
Thus, a somewhat different scenario explaining the arrival 
of fresh gas at the footprints of the GMLs may be required 
(Riquelme et al. 2010). Such may be the idea (Morris 2006; 
Torii et al. 2010) that the loops are magnetically-floated but 
accrete gas by driving shocks into the halo above the GC in 
which Hi reaches sufficient densities to efficiently cool and 
condense into H2- Of course, this general idea might be ex- 
tended to cover any sort of outflow that drives a sufficiently 
strong shock. 

Hi self-catalysed: Finally, we note that the Kelvin- 
Helmholtz instability will strip gas off the (relatively) cool 
and metal-rich Hi clouds ejected into the halo. Subsequent 
mixing of this stripped gas in the turbulent cloud wakes with 
in situ plasma allows for the cooling and condensation of 
the latter (allowing a phase transition from plasma to Hi) 
and the subsequent accumulation of some amount of this 
less astrated material (Marinacci et al. 2010, 2011; Binney & 
Fraternali 2011). Thus, not only does the fountaining of the 
cold, entrained gas present a mechanism for the net accretion 
of Hi in each cycle, it also potentially explains how new, rel- 
atively nuclear unprocessed matter finds its way into the GC 
system. Following the modelling of Marinacci et al. (2011) 
of this process in the disk of the Galaxy, net accumulation 
of gas of ~ 20% per fountain cycle is suggested; we thus 
plot 0.2 X Mcoid in fig. 21 as a rough estimator of the mass 
flow in this accretion channel. At the back-of-the-envelope 
level, this channel can supply a significant fraction - in fact 
saturate for high Mjn in the Esn[Mzams] scenario - the 
net accretion rate demanded by star-formation and the hot 
outflow, Mff^^^; clearly modelling particularised to condi- 
tions typical for the GC/inner Galaxy are required here to 
render this conclusion confident, however. 

Any of the three channels proposed above for accretion 
of matter out of the halo is of potential interest from a num- 
ber of points of view: i) as stressed, these mechanisms allow 
for the admixture of relatively pristine gas as demanded by 
observations - in general, this gas acts to counter-balance the 
effect of the high levels of astration in the GC region, thereby 
keeping overall metallicity somewhat in check; ii) such niech- 

We thank Richard Tuffs for raising this possibility. 



anisms show, in general, how the star-formation processes 
might be self-sustaining in the sense that the star-formation- 
driven outflow actually catalyses the further accretion of 
gas; and iii) following from this point, such a mechanism 
may explain how the long-term stability of the system - in 
particular, a star-formation rate apparently held rather con- 
stant over many dynamical times - is achieved (cf. Binney 
& Fraternali 2011; Marinacci et al. 2011). 

5.2 Strong, in situ magnetic field 

As remarked, our modelling shows that the in situ magnetic 
fleld is large ~ 100 /^G, potentially large enough to have 
important dynamical effects in the region. How is such a field 
established and - in light of the losses due to advection and 
turbulent diffusion (Beck et al. 1999) - maintained? Here a 
compelling mechanism is the long timescale accumulation of 
magnetic field lines frozen into gas accreted on to the region 
along the disk of the Galaxy (Chandran et al. 2000). This 
process produces a large-scale magnetic field orientation that, 
in agreement with observations, is predominantly vertical: 
compression amplifies the disk-perpendicular field component 
and advection of frozen-in field in the outflow (such as we 
have modelled) and/or ambipolar diffusion removes plane- 
parallel field components. Chandran ct al. (2000) estimate 
that the GC should accrete a total mass of ~ 3 x 10^ Mq over 
its lifetime (dependent on the amplitude of the assumed pre- 
Galactic field) in order to explain its present-day magnetic 
phenomenology; fig. 23 shows that this requirement is nicely 
matched for the regions of parameter space we favour (though 
we note that Chandran et al. 2000, favour a somewhat bigger, 
large-scale field amplitude than pointed to by our analysis). 

5.3 Implications of expelled, frozen-in magnetic 
field 

Figs. 7 and 8 point to magnetic energy losses at a rate of 
Xl, 3 X 10^^ erg/s. If sustained over 10 Gyr, these would 
imply the injection of ~ 10^^ erg by the outflow into 
the Fermi Bubbles, corresponding to a magnetic fleld am- 
plitude of 40 ^.G over the entire volume of the Bubbles; 
this is, however, a severe overestimate of the likely final 
fleld amplitude once the system has relaxed via reconnec- 
tion and adiabatic losses are accounted for (e.g., Braith- 
waite 2010). Magnetic field reconnection generically pro- 
vides for the local injection of energy in the Bubbles, de- 
livering heat and perhaps even providing for the accelera- 
tion of a non-thermal particle population in a distributed 
fashion in the Bubbles^^. Following Braithwaite (2010), we 
may roughly calculate t he amplitude of the fina l, relaxed 
field, Bf, as Bf ~ Bout\/ {r-out/rBub) {pBub/pout) ~ 100 pG 
X ^(110 pc/4.5 kpc) (0.01 cm-VO.l cm-3) ~ 5 ^G (here 
'oiit' denotes outflow and 'Bub' is for Bubble). This is charac- 
teristic amplitude over the entire volume of the Bubbles: we 
expect the field closer to the plane to be stronger in general 

Though note that the maximum rate of reconnection energy 
injection is a factor of a few short of being able to maintain 
the temperature of the plasma injected into the Bubbles: these 
radiative losses are sustained in steady state by the thermal power 
of the plasma itself after it is injected hot at the GC 
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(consistent with the fact that the non-thermal WMAP haze 
has now been shown to fairly sharply cut-off for b ^ 35°: 
Dobler 2011) and there to be local regions of enhanced field 
as demanded by the microwave polarisation observations 
(Jones et al. 2012). We note that within our model (Crocker 
& Aharonian 2011) that the Fermi Bubbles' 7-ray emission 
is supplied by hadronic processes, we have an expectation for 
the luminosity of the Bubbles due to synchrotron emission 
from secondary electrons and positrons (created in the same 
pp collisions delivering the 7-rays); to supply the observed 
non-thermal, microwave emission detected as the WMAP 
haze a ~ 10 fiG field is required. A further important effect 
of the injected magnetic field will be to stabilise the Bub- 
ble surfaces against fiuid instabilities (e.g. Gourgouliatos & 
Lyutikov 2012). 

5.4 Total masses ejected 

Fig. 23 also demonstrates something rather interesting about 
the total mass of hot gas processed through the system: 
this is (1 - 3) X 10^ Mq over the assumed 10 years for 
the favoured range of Min- In principle, there are three 
possible fates for this gas: i) it is ejected to infinity; ii) it 
rains back on to plane away from the GC; or iii) it falls 
back on the GC itself. Here ejection to infinity (in a true 
wind) can be firmly rejected as already noted. The exact 
proportion of material that falls back on to the GC will 
depend on the details of the interaction of the outflowing 
plasma with the differentially-rotating bulge (Marinacci et al. 
2011); we expect a non-negligible fraction of this low-angular 
momentum material to do so (see below). 

5.5 Relation to Fermi Bubbles 

Evident from fig. 26 is that the power going into the locally- 
accelerated cosmic ray proton particle population is almost 
invariant at ~ 10'^^ erg/s across the region of parameter space 
that fits the non-thermal data well. We nowhere constrain our 
modelling to reproduce this result - it is a prediction. This 
power is precisely enough to sustain the ~GeV 7-ray emission 
from the Fermi Bubbles assuming a saturation^'^ situation 
(Crocker & Aharonian 2011). It is also sufficient to supply 
the total few x 10^® erg enthalpy of the Bubbles were this 
power delivered for the multi-Gyr timescale independently 

Here 'saturation' implies that i) the system has 
reached steady state and that ii) the loss timescale of pp 
collisions is the smallest relevant timescale in the prob- 
lem, i.e., the Bubbles constitute thick targets to the in- 
jected cosmic ray protons because the pp loss time is, in 
particular, shorter than the escape time from the Bub- 
bles for protons of the relevant energy to generate ~GeV 
7-rays. If these conditions apply and cosmic rays remain 
co-entrained with the parcels of gas and magnetic field 
in which they are injected into the base of the Bubbles 
by the outflow, the volumetric emission from the Bub- 
bles becomes independent of the gas density in each such 
parcel (unless this is so low that the pp loss time be- 
comes longer than the age of the Bubbles). As detailed 
by Crocker & Aharonian (2011) this effect implies that 
the overall 7-ray emission from the Bubbles is rather uni- 
form as required by the observations. 



suggested by the long pp loss time (and the requirement 
that the steady-state be reached, achieving saturation; note 
that contributions from the outflowing plasma and magnetic 
fleld to the total enthalpy of the Bubbles would further 
decrease this timescale). This is surely not a coincidence. 
Also signiflcant is the fact that we explicitly conflrm (see 
flg. 7) that the outflow injects thermal plasma power into 
the base of the Bubbles easily sufflcient to explain their 
thermal emission of ~ 2 X 10^^ erg/s into the 0.5-2 keV 
energy band (Almy et al. 2000; Snowden et. al. 1997) or an 
inferred (temperature-dependent) bolometric luminosity of 
~ 10"° erg/s (Crocker & Aharonian 2011). 

The hot mass outflow suggested by our modelling (§5.4) 
is easily sufflcient to supply the ^ 10* Mq of plasma in the 
Fermi Bubbles over 10 Gyr but there is some sensitivity to 
the particularities of our modelling in this conclusion because 
the hot outflow is dominated by the centrally-loaded and 
heated ISM material. We can render the conclusion that GC 
star-formation can supply the plasma required to flU-out 
the Bubbles robust, however, by considering only the mass 
flux in the 'directly'-injected plasma (composed of supernova 
ejecta and stellar wind matter and neglecting the entrained 
ISM gas: there is little play in this number). Fig. 27 shows 
directly-injected plasma mass flux; even it alone is sufficient 
to sustain the mass growth rate of the Fermi Bubbles we 
have previously inferred (Crocker & Aharonian 2011) within 
our hadronic model (conclusion is robust with respect to 
both Min and choice of Esn and really just depends on the 
average SFR and the time over which it has been sustained), 
i.e., to supply ~ 10* Mq of gas over ~10 Gyr. 

Let us now proceed under the assumption that the GC 
outflow is feeding the Bubbles and that our model's treat- 
ment of mass loading of the hot outflow is correct. Then, if 
Mg °*10^ Mq of gas is ejected in the GC's hot outflow over 
*io^Gyr 10 Gyr, and fretMl;'°*lO^ Mq of gas is returned to 
the GC (after each cycle of the hot material from GC, into 
the halo and back again) and given the mass of the Fermi 
Bubbles if Mi'^W^ Mq (Su et al. 2010; Crocker & Aharo- 
nian 2011), the (initially) hot gas must cycle through the GC 
10 {fretMl^°^/M['^) times over i"-'^*, implying a cycle time 
of 1 X [t["(l%,^Mi'^ /ifretM^'")] Gyr. Self-consistently, this 
timescale is very similar to the cooling time of the ~ lO'' K, 
~0.01 cm~^ plasma in the Bubbles. The implication of this 
reasoning is that the gas in the Bubbles is slowly churning 
and, while the structures themselves may persist for ~ 10 
Gyr, individual parcels of gas injected into Bubbles remain 
aloft for only a fraction of this time with, presumably, a 
trajectory involving relatively high-speed injection in the 
outflow, a period of neutral or positive buoyancy relatively 
high in the halo (Rodriguez-Gonzalez et al. 2009), and flnally, 
condensation into Hi and fall-back to the GC once sufficient 
cooling has taken place. Zech et al. (2008) have previously 
interpreted data on high velocity clouds at inner Galactic 
longitudes as representing exactly such a circulation. 

If the Bubbles really are held aloft by buoyancy they 
must have an interior matter density less than the exter- 
nal halo plasma, implying a density contrast satisfying 
rj = pext/ pBub > 1- On the other hand, the Bubbles' asymp- 
totic ascent speed is controlled by t) with a larger density 
contrast implying a faster ascent (at least until additional, 
dynamical effects entrain signiflcant amounts of material 
into the Bubbles; e.g.. Pope et al. 2010). Following Hinton et 
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Figure 26. Model predictions for the power being fed into freshly- 
accelerated protons for the case (blue, solid) of i?sjv = 10^^ erg; 
and (purple, dashed) Esn an increasing function of Mzams as 
described in S3. 5. 
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Figure 27. Hot mass flows due to 'directly' injected material (i.e., 
supernova ejecta and stellar winds) for the case (blue, solid) the 
case of Egjq = 10^^ erg; and (purple, dashed) E^n an increasing 
function of Mzams as described in §3.5. The blue band shows 
the mass growth rate inferred (Crocker &: Aharonian 2011) for the 
Fermi Bubbles. 



al. (2007), our modelling allows us to calculate the density 
contrast assuming adiabaticity as 

7—1 ^ Loutrup . . 

■q = , — (6) 

7 M kg Te^t 

where ^ ~ 0.61 (assuming for simplicity a plasma of solar 
composition), the temperature of the halo plasma is T^^t ~ 
(3 - 10) X 10** K (Almy et al. 2000), and Lout and M are 
the mechanical power and mass injection rate of the outflow, 
respectively. In fig. 28 we calculate r] for the case of Esn ~ 
10^^ erg, taking the mass growth rate of the Bubbles to 
be given by Mhot (therefore ignoring mass drop-out due 
to cooling and condensation). Note that the star-formation- 
driven outflow we model naturally predicts an initial r; ~ 
(1 — 10); this is in contrast to AGN-driven outflows where 
initial density contrasts ?7 ^ 1 are generically expected (e.g., 
Hinton et al. 2007). We further find that 77 — > 1 (from above) 
in the range of Min preferred by other indicators predicting 
that - without fine-tuning - we expect the buoyant ascent 
of the Bubbles to be slow. 



5.6 Model-derived star formation and supernovae 
parameters 

On the basis of our modelling we can calculate - at least 
at the population level - various parameters characterising 
the region's star-formation and resultant supernovae. Fig. 29 
shows that our modelling suggests that as the favoured Min 



-1.0 -0.8 -0.6 -0.4 -0.2 
Log[A/,N/M /yrl 



0.0 



0.2 



Figure 28. Density contrast between external, halo plasma and 
Bubble interiors as filled by the outflow we model. The blue band 
is for an assumed halo plasma temperature of 10^'^ K and the red 
for 10^ K with the upper and lower range to each band given by 
taking the outflow fluid to be, respectively, a non-relativistic fluid 
(7 = 5/3) and a relativistic fluid (7 = 4/3). In this calculation 
adiabatic evolution is assumed and the mass flux into the Bubbles 
is just as given by Mhot ignoring mass drop-out due to cooling 
and condensation. 
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Figure 29. The lower Mzams cut-off - i.e., the fragmentation 
mass — as a function of Min- Note that at the very top of the 
allowed Min ranges the modelled fragmentation mass is butting- 
up against the (purely) numerical constraint that this parameter 
be larger than ~ 0.3 M©. 



parameter range (for either Esn case) is approached from 
below we see a transition from the modelling preferring 
high (1.2 M0) to low (i.e., conventional: 0.07 Af©) values 
of the fragmentation mass (the lower integration limit of 
the Mzams parameter). This is an interesting coincidence 
but numerical experiments we have performed indicate that 
our modelling does not seem to confer particularly strong 
sensitivity to this parameter. 

Fig. 30 demonstrates that, coincidentally or not, the 
modelled stellar wind losses within our scenario - as informed 
by a parameterization of the results presented by Meynet 
& Maeder (2003) for single star evolution - match rather 
well with the total stellar wind losses of ~ 0.01 Mq from the 
region determined by observations (see Crocker et al. 2011b, 
and references therein). 
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Figure 30. Stellar mass loss predicted by our modelling employing 
the prescription for single-star stellar winds presented by Meynet 
& Maeder (2003) for the cases of (blue) i?sjv = 10^^ erg (constant 
with respect to Mzams) ^^'^ (red) Egjv an increasing function of 
MzAMS ^ described in §3.5. The upper and lower bounds of each 
band are for the case, respectively, of stars rotating at 300 km/s 
and not rotating. The dashed line shows our previous estimate 
(Crocker et al. 2011b) for all GC stellar wind losses, O.OlMQ/year. 
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Figure 31. The total molecular gas mass in the GC region inferred 
from the Kennicutt (1998)-Schmidt relation and the modelled SFR 
for the case of (solid, blue) Egi^ = erg (constant with respect 
to Mzams) ^i^d (dashed, purple) i?sjv increasing function of 
Mzams ^ described in §3.5. 



Given the calculated amount of star-formation going 
on in the region, we can calculate an expectation for the 
region's total gas mass on the basis of the Kennicutt (1998)- 
Schmidt relation (see fig. 31); this expectation matches rea- 
sonably well to the estimated gas mass of 3 X 10^ Mq of 
the region (e.g., Molinari et al. 2011). In fact, if we demand 
that the GC be typically-efficient at turning gas into stars. 
Mi jv 1 Mq /year and Mjjv ^ 3 Mq /year are implied for 
the cases of Esn = 10^^ and Esn[Mzams], respectively. 
This is another (though perhaps rather weak) piece of evi- 
dence supporting Min ^ 1 Mo/year. 



5. 6. 1 Supernovae 

Fig. 32 shows that the modelled supernova rate in the region 
is 0.08-0.16 per century (corresponding to an expected time 
between supernovae of 630-1300 years); this is at the upper 



Figure 32. Supernova rate for the case of (solid, blue) Egff = 10^^ 
erg (constant with respect to Mzams) ^^<i (dashed, purple) i?sjv 
an increasing function of Mzams described in §3.5 and (dotted, 
yellow) the 'hypernova' rate, i.e., the number rate at which stars 
of Mzams > 26.5 - for which our parameterization satisfies 
Esn [Mzams] — 10^^ erg - are formed. 



end of the range we previously determined (viz. 0.02-0.08 
per century; see Appendix B2.2 of Crocker et al. 2011b), 
but certainly consistent with gross estimates (e.g., that the 
GC is responsible for a ~few-10 % of the Galaxy's massive 
star formation) and not in significant conflict with an upper 
limit of ~0.1 /century obtained from counting point-like radio 
source pulsar candidates (Lazio & Cordes 2008; Deneva et al. 
2009; also see Wharton et al. 2011). Much of the difference 
between our current supernova rate determination and our 
previous estimates is due to our adoption of a floating, lower 
cut-off to the fragmentation mass which, as shown above, 
the modelling prefers larger than the conventionally-adopted 
~ 0.07 M0 over much of the Mim parameter space. 

Also shown in this figure is the rate at which stars more 
massive than 26.5 Mq are produced which - in our adopted 
parameterizaton of Esn[Mzams] - is equivalent to the rate 
of SN explosions releasing more than 10^^ erg mechanical 
energy, i.e., hypernovae. This is similar to the total rate of 
broad-lined SNIbc expected in the Galaxy (adopting 2 core 
collapse supernovae per century, Diehl et al. 2006, of which 
1-2 % are SNIbc BL; e.g., Smith et al. 2011). 

Our modelling also allows us to calculate the mass 
ejected by each supernovae into the ISM (given certain as- 
sumptions: see fig. 33). We can compare these modelled 
masses with supernova ejecta masses inferred from observa- 
tions of supernova fight curves: 

SNIbc: With reasonable assumptions of average photo- 
spheric velocities Drout et al. (2011) have estimated the 
mean ejecta masses for SN lb, Ic, and Ic (broad line) types 
from photometric observations of a total of 25 Ibc super- 
novae detected within 150 Mpc with the robotic Palomar 
60-inch telescope. Consistent with the picture that either 
stellar winds, massive eruptions, or binary mass transfer (or 
some combination of these; e.g.. Smith et al. 2011) removes 
most of the envelopes of the massive progenitors of such 
explosions, Drout et al. (2011) determine mean ejecta masses 
for SN lb and Ic supernovae of ~ 2 A'Iq (and ~ 5 AIq for 
the more energetic, broad-line explosions); this result should 
be relatively unbiased by selection effects towards unusually 
bright SN Ibc^'*. 

SNII: Unfortunately, we are not aware of such an unbiased 
mean ejecta mass determination for SNII supernovae being 



Maria Drout, private communication 
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Figure 33. Mean ejecta mass from GC core-collapse supernovae 
assuming an invariant 1.4 Mq remnant produced by each SN 
and for different assumptions concerning annual mass loss rate 
for GC stars (which could, in principle be any or all of stellar 
winds, binary mass transfer, and eruptive mass loss) as indicated 
by the key with 'M&M' denoting the IMF-integrated effect of the 
single star stellar winds inferred from the prescription presented 
by Mcynct & Macdcr (2003) for non-rotating and stars initially 
rotating at 300 km/s, respectively. 
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Figure 34. Inferred rate of SNIb and SNIc explosions as a fraction 
of all core-collapse supernovae in the GC with different assump- 
tions for the total mass-loss rate experienced by all GC stars (i.e., 
through all channels, potentially including stellar winds, binary 
mass transfer, and eruptive mass losses). The Ibc/CCs„ith 2011 
horizontal line shows the cosmological average 26% contribution 
the sum of all SNIb and Ic make to the total core-collapse rate as 
inferred by Smith ct al. (2011). 



presented in the literature. With this in mind, from fig. S3 
from the Supplementary Information to the work by Perets 
et al. (2010), we find a rough mean SN II ejecta mass of 16 

Mq; from the data presented by Utrobin & Chugai (2011) 
on SN IIP, we find a mean 19.7 Mq ejecta mass. These latter 
determinations are quite consistent with SN II pregenitor's 
experiencing little mass loss at all during their lifetimes: with 
our assumed Kroupa IMF, the mean mass of stars formed 
between 8 and 150 Mq is 20.7 Mq, neglecting the mass of 
compact remnant, or 17.4 Mq, if a 1.4 Mq neutron stars 
is formed by all progenitors with Mzams between 8 and 20 

Mq and a ~ 8 Mq black hole is formed for more massive 
progenitors (Belczynski et al. 2011). Of course, these averages 
are presumably affected by selection biases as noted. 

Broadly, we find ejecta masses consistent with some ad- 
mixture of different SN types (fig. 33). In principle, we can 
actually infer the relative fractions of SNII and SNIbc types 
in the GC environment on the basis of our modelling and 
adopting the observational ejecta mass averages given above: 
see fig. 34 which shows the minimum fraction of SNIb and 
SNIc supernovae with respect to all core-collapse supernovae 
in the GC environment given the ejecta mass determined 
in our modelling and assuming that the progenitors of SNII 
experience negligible mass loss and form 1.4 Mq compact 
remnants. Adopting our fiducial number for the mass loss due 
to all stellar winds in the region, a SNIbc fraction ^ 50% is 
indicated, substantially larger than the cosmological average 
of 26% for all SNIb and Ic's arrived at by Smith et al. (2011). 
Allowing for other sources of mass loss (binary mass transfer, 
eruptive mass loss) would indicate an even higher fraction. 
These determinations are susceptible to many uncertainties, 
however, and we leave detailed treatment of this aspect of 
the problem to future work. For the moment we remark it 
is interesting that our modelling implies that the fraction of 
SNIb and SNIc events in the GC is higher than the cosmo- 
logical average. This is consistent with the finding that this 
fraction tends to increase in higher metallicity environments 
(e.g., Boissier & Prantzos 2009; Modjaz et al. 2011) and/or 
towards the brightest star-forming regions of gala^des (e.g., 
Leloudas et al. 2010). 
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Figure 35. Total numbers of compact stellar remnants formed 
over an assumed 10 Gyrs in the GC. Cases are given in the key; 'NS' 
denotes neutron stars and 'BH' denotes black holes. We assume 
that all stars 8 < Mzam S I ^0 < 20 form neutron stars and all 
stars Mzams > 2OAf0 form black holes. 



5.6.2 Stellar remnants 

Our modelling gives us a rough handle on the total number of 
compact remnants produced by GC star formation. Inspired 
by the findings of Belczynski et al. (2011), we assume that 
stars with 8 < Mzams / Mq < 20 form neutron stars and all 
stars Mzams > 20 Mq form black holes (with population- 
mean mass ~ 8 Mq ) . These assumptions would indicate 
(fig. 35) the formation of ~ lO'^ neutron stars and ~ 4 x 10® 
stellar mass black holes over the assumed 10 Gyr age of the 
system (with more total mass in black holes) 



Some fraction of the formed neutron star population will be 
both retained in the region and 'recycled' into the /is period range 
by spin-up due to accretion from a binary companion. In turn, 
some fraction - depending on beaming and other effects - of this 
milli-second pulsar (MSP) population can be expected to radiate 
at ~GeV energies (Abdo et al. 2010; Malyshev ct al. 2010). As we 
have previously remarked (Crocker ct al. 2011b), the ~ GeV flux 
data covering the region has a spectrum morphologically consistent 
with being dominated by such MSPs (which would have to amount 
to a few hundred in number). Given the many large uncertainties 
(particularly surrounding the recycling fraction), it is currently 
difficult to convert the neutron star formation rate (and assumed 
timescale) arrived at in our modelling to an accurate prediction 
of the number of extant and detectable MSPs. The analysis of 
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5.7 Implifications for picture of GC as 
starburst-like environment 

In disk environments - in botli our own Galaxy and in exter- 
nal, 'normal' spiral galaxies - it is generally held that ~ 10% 
of the initial kinetic energy delivered by a SN explosion is 
thermalized (Strickland & Heckman 2009, and references 
therein). In the GC environment we find, in contrast, a situ- 
ation much more akin - from this perspective - to starburst 
environments (Heckman et al. 1990). For instance, the recent 
detailed analysis of the starburst nucleus of M82 by Strick- 
land & Heckman (2009) points to a fraction approaching 
unity of each SN's energy being lost into heating or moving 
the ISM. Very hot plasma temperatures approaching 10** K 
also seem to be prevalent in such systems. 

As we have seen, for the GC, something similar seems to 
be happening. Prima facie, given the large amount of ambient 
molecular gas and the consequently large volumetric average 
gas densities in these systems, this appears surprising: why 
isn't more energy lost radiatively (cf. fig. 7)? 

Certainly, the volume filling factor of the densest cores 
of this gas (where most of the mass is tied up), is very small 
(e.g.. Oka et al. 2005). Moreover, the localisation of GC star- 
formation and its temporally-extended, continuous nature 
imply that the conditions in the GC resemble a sort of steady- 
state super bubble - or, indeed, starburst - where previous 
generations of stars and/or supernovae have presumably 
riven the ambient dense gas with cavities, blow-outs and 
channels (cf. Westmoquette et al. 2009; Contini 2011). These 
conditions mean shocks sent out into the ISM have a reduced 
opportunity to interact with dense ISM phases (that would 
then efficiently radiate) and, in fact, the energy injected by 
supernovae is efficiently thermalized as demanded (cf Higdon 
& Lingenfelter 2005). 

The same conditions mean that GC supernova remnants 
apparently remain good cosmic ray accelerators (assuming 
the EsN = 10®^ erg prescription is correct). The cosmic rays, 
moreover, only ever 'see' a small fraction of the total ambient 
gas mass (Crocker et al. 2011b); both plasma and cosmic 
rays are apparently blown out of the system so quickly that 
they they do not penetrate into the dense gas cores where 
stars are being formed. 



6 DISCUSSION 

6.1 Is the 'very hot' plasma real or not? 

We have already seen (fig. 25) that the Esn = 10^^ erg case 
indicates a plasma temperature ^ 10^ K over the favoured 
MiN parameter space; on the other hand, the Esn[Mzams] 
case provides for a temperature approaching ~ 7 x lO'^ K 
at the lower end of its favoured AIjn range. Moreover, if 
we constrain our modelling so that the modelled plasma 
temperature is ^ 5 x 10^ K we find good fits to the data for 
even the Esn = 10^^ erg case, though, of course, the plasma 
temperature is no longer then a prediction of the modelling. 
It would seem yet premature, therefore, to rule the idea of 
the very hot plasma out on the basis of our modelling. 

Wharton et al. (2011) seems to broadly show that our numbers 
are reasonable given the play in the relevant parameters. 



Two other facts should also be kept in mind here. Firstly, 
we are modelling the plasma's temperature as it leaves the 
region of interest; this is after the centralised mass loading 
(i.e., the heating of swept-up, ambient ISM gas to plasma 
temperatures), also pointed to by our analysis, has taken 
place. This mass loading process has to take place over some 
physical scale (and may, indeed, continue once the outflowing 
gas has left the region). The natural temperature scale for 
the plasmas in the centres of star-forming regions before any 
mass loading has taken place (and assuming themalization 
efficiency /3 — >■ 1, not too much at variance with what we 
find) is ~ 3 X 10* K (e.g., Veilleux et al. 2005). So, in 
principle, the X-ray observations may be revealing the very 
hot, central temperature of the plasma before too much mass 
loading takes place (as the gas is driven out of the system) . 
Indeed, emission from the putative very hot component is 
concentrated on somewhat smaller size scales than the full 
200 X 80 pc^ region modelled here (Yamauchi et al. 1990; 
Belmont et al. 2005). 

Previous work on the question of the existence 
of the putative very hot plasma (e.g., Muno et al. 
2004) has focused on the energetic difficulties ap- 
parently encountered were the plasma real and es- 
caping at approximately the sound speed because 
of the inadequacy of known power sources to sus- 
tain the thermal energy carried off by the plasma. 
Setting aside the fact that our modelling apparently 
favours a lower temperature for the escaping diffuse 
plasma anyway, it is worth emphasising here that 
i) the mass loading and entrainment processes we 
have modelled here slow down the outflow signifi- 
cantly with respect to the sound speed and ii) the 
supernova rate we find — close to 1 per 1000 years — 
is by construction sufficient to sustain the energy losses 
represented by the outfiow. 

Note that other pieces of evidence continue to suggest 
the plasma is real: the suggestive pressure equilibrium that 
would exist between such a plasma and the other GC ISM 
phases (Spergel & Blitz 1992); the fact that such very hot, 
diffuse plasmas are seen in the centres of active star-forming 
systems like M82 that are driving 'super' outflows (Strickland 
& Heckman 2009) and, strengthening this connection, the 
fact that the 6.7 keV Fe Hea line flux measured from the 
GC (~ 4 X 10^* erg/s as inferred from table 2 of Koyama 
et al. 2007) is exactly as expected (see §5.5 of Strickland & 
Heckman 2009) given the GC SFR and overall size of the 
star-forming region. 

On the other hand, as explained above, we find evidence 
for (relatively) cool, entrained gas in our outflow, apparently 
with rather large flUing factor, consistent with the envelope 
H2 phase identifled by and other measurements; it has 
been claimed (Goto et al. 2008) that it is difficult to reconcile 
the existence of this phase with the existence of a (necessarily) 
large filling factor very hot plasma. In summary, it seems 
that for the moment the jury remains out on the important 
question of whether the very hot, diffuse plasma exists or 
not. 
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6.2 Is the system really in steady state? What 
about the contribution of individual 
star-bursts? 

The steady state modelling presented here indicates that 
~ lO'^ '^ core-collapse supernovae have occurred in the GC in 
the post-infant Galaxy. This number is ultimately determined 
by the secular accretion of gas on to the region which has 
led to the formation of ~ 10^ Mq of stars in the nuclear 
bulge. Of course there are stochastic variations in the GC 
star-formation rate. According to the analysis of Sjouwerman 
et al. (1998) of their observations of oxygen-rich, cool giant 
stars, the most recent, significant burst of SF in the GC 
occurred approximately 1 Gyr ago and should have resulted 
(Higdon et al. 2009) in ~ 10^ core-collapse supernovae. Thus 
even this relatively dramatic event in the GC's history only 
contributes at the ^ 10% to the total star-formation (and 
resultant energy-deposition) history of the region. 

Observations of today's GC super stellar clusters of 
stochastic variation in the star-formation rate on much more 
recent timescales, few x 10^ Myr, but the total mass of stars 
formed in these systems/events, ^ 10^ Mq, is much less than 
the ~ 10^ Mq of stars represented by the nuclear bulge. Such 
events will only be associated with the ejection of significant 
amounts of energy if they are accompanied by substantial 
jumps in the accretion rate on to Sgr A* (cf. Zubovas et al. 
2011). 

6.3 How is the long-term stability of the SFR 
maintained? 

The long-term star formation scenario requires that the GC 
continues to accumulate gas over a similar timeframe; prima 
facie this requires that a driving stellar bar remains extant, 
too, over the requisite multi Gyr timescale (Serabyn & Morris 
1996). The fact that bars are common in spiral galaxies 
(Eskridge et al. 2000) is consistent with this requirement. 
On the other hand, we have listed a number of mechanisms 
that may allow the GC to pull gas directly out of the halo. 
The operation of some such mechanism seems to be required 
by the detection of gas that has undergone relatively little 
astration in the GC region and may imply a 'guaranteed' 
accretion rate that means that GC is not completely at the 
mercy of conditions prevailing on much wider scales in the 
disk to sustain its star-formation. 

The scenario investigated here of a mass-loaded, sub- 
sonic outflow is interestingly reminiscent of the situation 
identified by Tang ct al. (2009) who model long-lasting feed- 
back due to outflows from galactic nuclei (driven by Type 
la supernovae) wherein there is little active star-formation. 
These authors find that, by dropping out cool, entrained 
gas at large radii, such outfiows are stabilised to the extent 
that they remain essentially stationary for a few Gyr. As 
we have previously argued (Crocker & Aharonian 2011) and 
continue to find in the research described here, the evidence 
for the GC system and the outflow it drives is that a similar 
situation of quasi-stationarity has been reached. 

Finally, Kim et al. (2011) in their numerical treatment of 
GC star-formation note that a negative feedback mechanism 
may be implied by the fact that GC star-formation peaks 
sharply in a region consistent in extent with the expected size 
of the X2 orbits (i.e., a smaller region than the whole CMZ) 



but seems to be relatively suppressed beyond this region 
(in the XI orbit region) despite the presence of gas. They 
postulate that heating of the relatively more diffuse gas in a 
UV photon field of ~ 30 eV cm~'^ - not unreasonable for the 
inner GC environment - may provide for this suppression. 
Remembering our finding above and previously (Crocker et 
al. 2011a, b) that cosmic rays do not seem to penetrate into 
the dense gas in the time they remain in the system, we 
note here that cosmic ray heating may equally-well affect 
(or at least supplement) the requisite negative feedback (cf. 
Suchkov et al. 1996; Yusef-Zadeh et al. 2007; Crocker et al. 
2011b, and references therein). Consistent with this, we find 
steady-state cosmic ray energy densities in our modelling of 
> 20 -40 eV cm"^. 



7 CONCLUSIONS 

We have modelled the mass and energy flows through the 
Galactic centre in a one-zone model and shown that star- 
formation - and resultant supernovae - in the GC are well- 
sufficient to drive the gross dynamics of the system and 
to explain its non-thermal phenomenology. None of this is 
to deny the importance of, e.g., stellar-radiation driving of 
gas dynamics in particular regions of the GC, but radia- 
tion pressure is not required to explain the large-scale mass 
movements we infer here. 

Current evidence does not seem strong enough for us 
to promote the empirically-inspired scaling of supernova 
mechanical energy with zero-age, main sequence progeni- 
tor mass that we have tested above; if anything, the 'stan- 
dard' EsN = lO^'^ erg assumption seems to work better in 
a number of cases (e.g., in suggesting a mass loading fac- 
tor closer to expectation for star-burst-like environments; a 
SFR closer to other, independent estimates; supplying the 
total stellar mass of the nuclear stellar bulge for less ex- 
treme values of Min)- In any case, as emphasised above, 
the EsN scaling we adopted certainly constitutes an upper 
limit to the true, population-averaged evolution of Esn with 
stellar mass. Likely, further modelling, probably involving 
other constraints, will be needed before a better handle on 
Esn[Mzams] can be arrived at. Given these considerations, 
we now particularise our discussion to the case of Esn = 10^^ 
erg. 

A number of indicators come together to suggest that 
our control parameter, M/jv - the total mass being fed into 
the system - has a lower limit at around 0.4 /year. We 
find a 2(7 upper limit on Mjn (for the Esn = 10^^ case) 
at 1.8 M0/year. We emphasise, however, that Min is an 
overestimate of the total mass flux accreting out of the plane 
of the Galaxy on to the GC. Our modelling suggests that 
there is an outflow of plasma and cosmic rays from the system 
that entrains cool gas. This entrained mass constitutes most 
of the mass flux but will fountain back on the GC. This 
naturally accounts for the recently-observed halo of warm 
molecular hydrogen found to be surrounding the CMZ and 
may self-catalyse the accretion of relatively pristine corona 
gas into the system. 

We flnd that our modelling robustly predicts an almost 
invariant 10'^^ erg/s for the power going into the freshly 
accelerated cosmic ray proton population in the GC region. 
We nowhere constrain our model to produce this result: 
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it emerges from the numerical minimisation procedure. As 
we have previously emphasised (Crocker & Aharonian 2011; 
Crocker et al. 2011b) this power is precisely enough to sustain 
the 7-ray emission from the Fermi Bubbles in a hadronic 
saturation scenario and allow for the inflation of the Bubbles 
against the pressure of the external medium in a few Gyr 
(probably assisted by the injected magnetic field and injected 
plasma) . Equally, we find that the modelling robustly predicts 
that the GC system injects ~ 10''* erg/s into hard-spectrum 
cosmic ray electrons; this is sufficient (Crocker et al. 2010a, 
2011b) to explain the non-thermal synchrotron radiation 
detected from the GC lobe and wider diffuse, non-thermal 
source region detected around the GC (LaRosa et al. 2005; 
Crocker et al. 2010a). 

Together with the evidence that the GC SFR has been 
quasi-stationary for Gyr timescales and that the outflow 
from the GC advects most non-thermal particles out of the 
acceleration region before they lose much energy, we consider 
our finding that the GC accelerates 10^^ erg/s in cosmic rays 
as a very strong indication that GC star-formation essentially 
explains the Fermi Bubbles. This argument is strengthened 
by the facts that the same star-formation processes can inject 
the plasma mass and thermal power required to fill-out the 
Bubbles and sustain their X-ray emission and to inject the 
magnetic fields that can stabilise the Bubble surfaces against 
fiuid instabilities, trap their cosmic ray and plasma contents 
for long timescales and explain their microwave synchrotron 
emission. 



The Galactic centre is not particularly distinguished in 
the night sky - but this belies its true activity: the many 
orders of magnitude of visual extinction arising from the 
column of dust we view it through is the reason for this 
extrinsic dimness. Equally, the GC is not particularly im- 
pressive as a non-thermal radiation source; its ~ GeV 7-ray 
luminosity, at few xlO^® erg/s, is an order of magnitude 
short of the 5-10% of the Galaxy's output one might guess 
on the basis of the amount of massive star-formation hap- 
pening in the system (the Galactic ~ GeV 7-ray luminosity 
is about 3 x 10^* erg/s as inferred from fig. 1 of Strong et 
al. 2010). As we have emphasised previously (Crocker et al. 
2011a; Crocker & Aharonian 2011; Crocker ct al. 2011b) this 
is because the cosmic rays accelerated in the region are mostly 
leaving before they can radiate. The radiation that these par- 
ticles finally do emit is writ large in the Fermi Bubbles: it is 
from these structures that we detect the ~ 10% of Galactic 
7-ray luminosity (i.e., 4 x 10^^ erg/s) that we expect on the 
basis of the GC's share of Galactic star-formation. 

The fact that we can detect the Fermi Bubbles at all is 
testament to the long-term stability of the GC as star-forming 
system. We have seen hints above as to how this stability can 
be established, in particular, how GC star-formation activity 
can be insulated from the vicissitudes of conditions in the 
Galactic disk: it seems that a minimal level of accretion 
in the GC system is self-catalysed. This is consistent with 
the presence of relatively pristine gas accreted out of the 
halo via a mechanism or mechanisms directly related to the 
star-formation-driven outflow (i.e., injection of cool, high 
metallicity gas and/or dust and/or shocks into the halo 
plasma) . 



We find complete consistency between the long-timescale- 
averaged power required to drive the GC system and the 
power injected by the star-formation we can infer is currently 
taking place in the system. We need not invoke periods of very 
bright AGN-type activity of the supermassive black hole to 
explain the dynamics of the GC. In fact, as a final speculation: 
these studies hint that the importance of sustained, GC star- 
formation is that it effectively erects a curtain wall around 
the SMBH, either using-up gas directly by creating new stars 
or blowing it away before much of it can reach Sgr A*. This 
prevents mass accretion at rates that would allow the SMBH 
to undergo phases of activity sufficient to heat the Galaxy's 
coronal gas to such high temperatures that further accretion 
on the Galactic disk would become impossible^*'. This, in 
turn, enables the long-term sustenance of disk star formation 
(cf. Binncy & Fratornali 2011). Such a mechanism would 
explain the emerging finding (Erwin & Gadotti 2011) - to 
which the Milky Way adheres - that the mass of a nuclear 
star cluster correlates with the total stellar mass of its host 
galaxy rather than the galaxy's bulge. 
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APPENDIX A: METALLICITY OF GC 
ENVIRONMENT 

There remain considerable uncertainties surrounding the 
metallicity of the GC environment with different observa- 
tional probes suggesting different results. Broadly, however, 
these different approaches suggest that metallicity of the envi- 
ronment is between solar and twice solar: Zq ^ Zee ~ 2^0. 

In more detail, from measurements of the nitrogen sur- 
face abundance of WN stars in the Arches cluster Najarro 
et al. (2004) report a solar metallicity. Grafener et al. (2011) 
have, however, claimed that these measurements actually 
represent a lower limit to the initial C + N + O abundance 



Figure Bl. An example modelled broadband spectrum for the 
region (showing non- thermal components only) , in this case for the 
control parameter M/jv = 1 Mq /year. The curves are as follows: i) 
solid blue: total synchrotron; ii) dashed blue: synchrotron from pri- 
mary electrons (almost indistinguishable from total synchrotron in 
this case); iii) dotted blue: synchrotron from secondary electrons; 
iv) solid purple: sum of 7-ray fluxes from all processes; v) dot-dash 
brown: pp — >■ neutral mesons — > 7-rays; vi) dashed green: inverse 
Compton by primary electrons; vii) dotted green: : inverse Comp- 
ton by secondary electrons; viii) dashed red: bremsstrahlung by 
primary electrons; and ix) dotted red: bremsstrahlung by secondary 
electrons. Note the ~GeV data points and the X-ray flux datum 
only constitute weak upper limits to the diffuse, non-thermal emis- 
sion (the X-ray datum is explained via thermal bremsstrahlung 
emission by the plasma, an emission process we also model). 

with a best-fit value around 2Zq. Recent spectral analy- 
sis of luminous cool stars (Cunha et al. 2007) and LBVs 
(Najarro et al. 2009) indicates an approximately solar Fe 
abundance and enhanced abundances of a-elements. Lier- 
mann et al. (2010), citing the Martins et al. (2008) survey 
of the Arches cluster, determine a slight metallicity enrich- 
ment: Zac = (1.3 - 1A)Zq. Mauerhan et al. (2010) in their 
Paschen-Q survey of the inner ~90 pc determine a domi- 
nance of WNL and WCL Wolf-Rayet subtypes over WNE 
and WCE subtypes. This is consistent with a metallicity for 
the region that is at least solar. 

From X-ray observations Wang et al. (2006) report a 
plasma Fe abundance I.SIq ^ solar; Koyama et al. (2007) 
and Nobukawa et al. (2010) report an upper limit at twice 
solar for the inner ~ 200 pc region. Borkowski et al. (2010) 
report super-solar metal abundances in the GC supernova 
remnant Gl. 9-1-0. 3 on the basis of X-ray line observations. 
On the basis of high-latitude UV spectroscopy Zech et al. 
(2008) report the existence of super solar metallicity gas 
likely associated with a GC wind or fountain. 

Finally, studies of radio recombination line emission 
from Hll regions reveal a clear gradient in metallicity with 
Galactic radius (Shaver et al. 1983; Afflerbach et al. 1996) 
and predict that the warm plasma phase of the GC ISM 
should be rather cool; radio recombination line studies (Law 
et al. 2009) of the GC reveal exactly this. 



APPENDIX B: EXAMPLE BROADBAND FIT 

We show here an example broad-band spectral fit obtained 
by our -minimization procedure. 



